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This paper presents a comprehensive analysis of simple models useful to analyze the growth of 
nanostructures obtained by cluster deposition. After detailing the potential interest of nanostruc- 
tures, I extensively study the first stages of growth (the submonolayer regime) by kinetic Monte- 
Carlo simulations. These simulations are performed in a wide variety of experimental situations 
: complete condensation, growth with reevaporation, nucleation on defects, total or null cluster- 
cluster coalescence .... The main scope of the paper is to help experimentalists analyzing their 
data to deduce which of those processes are important and to quantify them. A software including 
all these simulation programs is available at no cost on request to the author. I carefully discuss 
experiments of growth from cluster beams and show how the mobility of the clusters on the surface 
can be measured : surprisingly high values are found. An important issue for future technological 
' applications of cluster deposition is the relation between the size of the incident clusters and the size 

of the islands obtained on the substrate. An approximate formula which gives the ratio of the two 
sizes as a function of the melting temperature of the material deposited is given. Finally, I study 
the atomic mechanisms which can explain the diffusion of the clusters on a substrate and the result 
of their mutual interaction (simple juxtaposition, partial or total coalescence . . . ). 
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I. INTRODUCTION 



Growth of new materials with tailored properties is one of the most active research directions for physicists. As 
pointed out by Silvan Schweber in his brilliant analysis of the evolution of physics after World War II : "An important 
transformation has taken place in physics : As had previously happened in chemistry, an ever larger fraction of 
the efforts in the field [are] being devoted to the study of novelty [creation of new structures, new objects and new 
phenomena] rather than to the elucidation of fundamental laws and interactions [...] Condensed matter physics has 
indeed become the study of systems that have never before existed." 

Among these new materials, those presenting a structure controlled down to the nanometer scale are being ex- 
tensively studied B]?]. There are different ways to build up nanostructured systems Q : atomic deposition ||, 
£f\ , mechanical milling fjlC j , chemical methods jTjJTTl , gas-aggregation techniques [^ JT^JT^ ] . . . Each of these techniques has 
its own advantages, but, as happens with atomic deposition techniques, the requisites of control (in terms of character- 
&\ I ization and flexibility) and efficiency (in terms of quantity of matter obtained per second) are generally incompatible. 
0^ ■ As a physicist wishing to understand the details of the processes involved in the building of these nanostructures, 
I will focus in this review on a carefully controlled method : low energy cluster deposition pJ. Clusters are large 
"molecules" containing typically from 10 to 2000 atoms, and have been studied for their specific physical properties 
(mostly due to their large surface to volume ratio) which are size dependent and different from both the atoms and the 
, bulk material [Q,|l4|-[F^] . By depositing preformed clusters on a substrate, one can build nanostructures of two types 
: in the submonolayer range, separated (and hopefully ordered) nanoislands, and for higher thicknesses, thin films or 
cluster assembled materials (CAM) . The main advantage of the cluster deposition technique is that one can carefully 
control the building block (i.e. the cluster) and characterize the growth mechanisms. By changing the size of the 
incident clusters one can change the growth mechanisms Jl8],[l9| and the characteristics of the materials. For example, 
• i-H . it has been shown that by changing the mean size of the incident carbon clusters, one can modify the properties of 
the carbon film, from graphitic to diamond- like p0[ . 

This review is organized as follows. First, I present briefly the interest of nanostructures, both in the domain of 
nanoislands arranged on a substrate and as nanostructured, continuous films. I also review the different strategies 
employed to deposit clusters on a substrate : by accelerating them or by achieving their soft-landing. The scope of 
this section is to convince the reader that cluster deposition is a promising technique for nanostructure growth in a 



variety of domains, and therefore deserves a careful study. In Section II], models for cluster deposition are introduced. 
These models can also be useful for atomic deposition in some simple cases, namely when aggregation is irreversible. 
The models are adapted here to the physics of cluster deposition. In this case, reevaporation from the substrate 
can be important (as opposed to the usual conditions of Molecular Beam Epitaxy), cluster-cluster aggregation is 
always irreversible (as opposed to the possibility of bond breaking for atoms pl[ ) and particle-particle coalescence is 
possible. After a brief presentation of Kinetic Monte-Carlo (KMC) simulations, I show how the submonolayer regime 



1 



can be studied in a wide variety of experimental situations : complete condensation, growth with reevaporation, 
nucleation on defects, formation of two and three dimensional islands . . . Since I want these models to be useful for 
experimentalists, Section [v| is entirely devoted to the presentation of a strategy on how to analyze experimental 
data and extract microscopic parameters such as diffusion and evaporation rates. I remind the reader that a simple 
software simulating all these situations is available at no cost on simple request to the author. Section ^ analyzes in 
detail several experiments of cluster deposition. These studies serve as examples of the recipes given in Section [v] to 
analyze the data and also to demonstrate that clusters can have surprisingly large mobilities (comparable to atomic 
mobilities) on some substrates. A first interpretation of these intriguing results at the atomic level is given in Section 
|VII| , where the kinetics of cluster-cluster coalescence is also studied. The main results of this Section are that high 
cluster mobilities can be achieved provided the cluster does not find an epitaxial arrangement on the substrate and 
that cluster-cluster coalescence can be much slower than predicted by macroscopic theories. 

A note on terminology : The structures formed on the surface by aggregation of the clusters are called islands. 
This is to avoid the possible confusion with the terms usually employed for atomic deposition where the clusters are 
the islands formed by aggregation of atoms on the surface. Here, the clusters are preformed in the gas phase before 
deposition. I use coverage for the actual portion of the surface covered by the islands and thickness for the total 
amount of matter deposited on the surface (see also Table I) . 



II. INTEREST AND BUILDING OF NANOSTRUCTURES 



Before turning to the heart of this paper - the growth of nanostructures by cluster deposition - I think it is 
appropriate to show why one wants to obtain nanostructures at all and how these can be prepared experimentally. 
Due to the technological impetus, a tremendous amount of both experimental and theoretical work has been carried 
out in this field, and it is impossible to summarize every aspect of it here. For a recent and rather thorough review, 
see Ref. || where the possible technological impact of nanostructures is also addressed. Actually, several journals 
are entirely devoted to this field [^2| . The reader is also referred to the enormous number of World Wide Web pages 
(about 6000 on nanostructures), especially those quoted in Ref. (^3|. A short summary of the industrial interest of 
nanostructures [Hi and introductory reviews on the interest of "Nanoscale and ultrafast devices" (2^] or "Optics of 
nanostructures" |26| have appeared recently. 

There are two distinct (though related) domains where nanostructures can be interesting for applications. The 
first stems from the desire of miniaturization of electronic devices. Specifically, one would like to grow organized 
nanometer size islands with specific electronic properties. As a consequence, an impressive quantity of deposition 
techniques have been developed to grow carefully controlled thin films and nanostructures from atomic deposition JB| . 
While most of these techniques are complex and keyed to specific applications, Molecular Beam Epitaxy (MBE) |27]] 
has received much attention from the physicists |^8| , mainly because of its (relative) simplicity. The second subfield 
is that of nanostructured materials as thin or thick films, which show (mechanical, catalytic, optical) properties 
different from their microcrystalline counterparts P,p|, pl] , p9pC| ] . 

I will now briefly review the two subfields since cluster deposition can be used to build both types of nanostructures. 
Moreover, some of the physical processes studied below (such as cluster-cluster coalescence) are of interest for both 
types of structure. 



A. Organized nanoislands 



There has been a growing interest for the fabrication of organized islands of nanometer dimensions. One of the 
reasons is the obvious advantage of miniaturizing the electronic devices both for device speed and density on a chip 
(for a simple and enjoyable introduction to the progressive miniaturization of electronics devices, see Ref. [f?T|). But 
it should be noted that at these scales, shrinking the size of the devices does also change their properties, owing to 
quantun confinement effects. Specifically, semiconductor islands smaller than the Bohr diameter of the bulk material 
(from several nm to several tens of nm) show interesting properties : as their size decreases, their effective bandgap 
increases. The possiblity to tailor the electronic properties of a given material by playing on its size has generated 
a high level of interest in the field of these quantum dots [j32[ . But quantum dots are not the only driving force for 
obtaining organized nanoislands. Isolated nanoparticles are also interesting as model catalysts (see Refs. |33|,|34| and 
Chapter 12 of Ref. ||). Clearly, using small particles increases the specific catalytic area for a given volume. More 
interesting, particles smaller than 4-5 nm in diameter might show specific catalytic properties, different from the bulk 
p4j , f35t , although the precise mechanisms are not always well identified (Chapter 12 of Ref. ^|). One possibility is 
the increase, for small particle sizes, of the proportion of low coordination atoms (corners, kinks) whose electronic 
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(and therefore catalytic) properties are expected to be different from bulk atoms. For even smaller particles (1-2 
nm), the interaction with the substrate can significantly alter their electronic properties p6| . Recently, there have 
been attempts at organizing the isolated islands to test the consequences on the catalytic properties |37j . Obtaining 
isolated clusters on a surface can also be interesting to study their properties. For example, Schaefer et al. |3|] have 
obtained isolated gold clusters onto a variety of substrates to investigate the elastic properties of single nanoparticles 
by Atomic Force Microscopy (AFM). 

Let me now briefly turn on to the the possible ways of obtaining such organized nanoislands. Deposition of atoms 
on carefully controlled substrates is the main technique used presently by physicists to try to obtain a periodic array 
of nanometer islands of well-defined sizes. A striking example J3!| of organized nanoislands is given in Fig. These 
triangular islands have been grown on the dislocation network formed by the second Ag atomic layer on Pt(lll). 
Beautiful as these triangles are, they have to be formed by nucleation and growth on the substrate, and therefore 
the process is highly dependent on the interaction of the adatoms with the substrate (energy barriers for diffusion, 
possibility of exchange of adatoms and substrate atoms, . . . ). This drastically limits the range of possible materials 
that can be grown by this method. However, the growth of strained islands by heteroepitaxy is under active study, 
since stress is a force which can lead to order, and even a tunable order, as observed for example in the system 
PbSejPb X - x Eu x Te @ (see Refs. @|| for further details on stress). 



In this review, I will focus on an alternative approach to form nanoislands on substrates : instead of growing them by 
atom-atom aggregation on the substrate, a process which dramatically depends on the idiosyncrasies of the substrate 
and its interaction with the deposited atoms, one can prepare the islands (as free clusters) before deposition and then 
deposit them. It should be noted that the cluster structure can be extensively characterized prior to deposition by 
several in-flight techniques such as time-of-flight spectrometry, photo-ionization or fragmentation p3[ . Moreover, the 
properties of these building blocks can be adjusted by changing their size, which also affects the growth mechanisms, 
and therefore the film morphology [|l8Ul9| . A clear example of the possibility to change the film morphology by 
varying only the mean cluster size has been given a few years ago by Fuchs et al. |lj] (Fig. ||) and this study has 
been completed recently by Brechignac's group for larger cluster sizes [ [l9|| . There are several additional interests for 
depositing clusters. First, these are grown in extreme nonequilibrium conditions, which allows to obtain metastable 
structures or alloys. It is true that neither islands grown on a substrate are generally in equilibrium, but the quenching 
rate is very high in a beam, and the method is more flexible since one avoids the effects of nucleation and growth on 
a specific substrate. For example, PdPt alloy clusters - which are known to have interesting catalytic properties - can 
be prepared with a precise composition (corresponding to the composition of the target rod, see below) and variable 
size and then deposited on a surface H]. The same is true for SiC clusters where one can modify the electronic 
properties of the famous Cqq clusters by introducing in a controlled way Si atoms before deposition |^5| . This allows 
to tune within a certain range the properties of the films by choosing the preparation conditions of the preformed 
clusters. It might also be anticipated that cluster nucleation is less sensitive to impurities than atomic nucleation. 
Atomic island growth can be dramatically affected by them, as exemplified by the celebrated case of the different 
morphologies of Pt islands grown on Pt(lll) which were actually the result of CO contamination at an incredibly 
low level : 10 -10 mbar [j47|. Instead, clusters, being larger entities, might interact less specifically with the substrate 
and its impurities. There is still no systematic way of organizing the clusters on a surface. One could try to pin them 
on selected sites such as defects or to encapsulate the clusters with organic molecules before deposition in order to 
obtain ordered arrays on a substrate Q . 

B. Nanostructured materials 

Although my main focus in this review is the understanding of the first stages of growth, it is worth pointing out 
the interest of thicker nanostructured films (for a recent review of this field, see Ref. |g]). It is known P,p|j29|] that the 
(magnetic, optical and mechanical) properties of these films can be intrinsically different from their macrocrystalline 
counterparts. The precise reasons for this are currently being investigated, but one can cite the presence of a signifi- 
cant fraction (more than 10 %) of atoms in configurations different from the bulk configuration, for example in grain 
boundaries It is reasonable to suppose that both dislocation generation and mobility may become significantly 
difficult in nanostructured films |5). For example, recent studies of the mechanical deformation properties of nanocrys- 
talline copper [[l9| have shown that high strain can be reached before the appearance of plastic deformation. A review 
of the effects of nanostructuration on the mechanical response of solids is given by Weertman and Averback in chapter 
13 of Ref. Another interesting property of these materials is that their crystalline order is intermediate between 
that of the amorphous materials (first neighbor order) and of crystalline materials (long range order). It is given by 
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the size of the crystalline cluster, which can be tuned. For example, for random magnetic materials, by varying the 
size of the clusters, and consequently of the ferromagnetic domain, one can study the models of amorphous magnetic 
solids @. 



C. How can one deposit clusters on surfaces? 



After detailing the potential interests of nanostructures, I now address the practical preparation methods by cluster 
deposition. Two main variants have been explored. Historically, the first idea has been to produce beams of accelerated 
(ionized) clusters and take advantage of the incident kinetic energy to enhance atomic mobility even at low substrate 
temperatures. This method does not lead in general to nanostructured materials, but to films similar to those obtained 
by atomic deposition, with sometimes better properties. A more recent approach is to deposit neutral clusters, with 
low energy to preserve their peculiar properties when they reach the surface. The limit between the two methods is 
roughly at a kinetic energy of 0.1 to 1 eV/atom. 



1. Accelerated clusters 



The Japanese group of Kyoto University was the first to explore the idea of depositing clusters with high kinetic 
energies (typically a few keV) to form thin films The basic idea of the Ionized Cluster Beam (ICB) technique is 
that the cluster breaks upon arrival and its kinetic energy is transferred to the adatoms which then have high lateral 
(i.e. parallel to the substrate) mobilities on the surface. This allows in principle to achieve epitaxy at low substrate 
temperatures, which is interesting to avoid diffusion at interfaces or other activated processes. Several examples of 
good epitaxy by ICB have been obtained by Kyoto's group : Al/Si |5^] which has a large mismatch and many other 
couples of metals and ceramics on various crystalline substrates such as Si(100), Si(lll) . . . Molecular Dynamics (MD) 
simulations have supported this idea of epitaxy by cluster spreading ||| . The reader is referred to Yamada's reviews 
pl[ for an exhaustive list of ICB applications, which also includes high energy density bombardment of surfaces to 
achieve sputter yields significantly higher than obtained from atomic bombardment |p4| . 

However, the physics behind these technological successes is not clear. In fact, the very presence of a significant 
fraction of large clusters in the beam seems dubious J55||36]]. There is some experimental evidence |5l[] offered by 
Kyoto's group to support the effective presence of a significant fraction of large clusters in the beam, but the evidence 
is not conclusive. In short, it is difficult to make a definite judgement about the ICB technique. There is no clear 
proof of the presence of clusters in the beam and the high energy of the incident particles renders difficult any attempt 
of modelling. Kyoto's group has clearly shown that ICB does lead to good quality films in many cases but it is not 
clear how systematic the improvement is when compared to atomic deposition techniques. 

Haberland's group in Freiburg has developed recently a different technique called Energetic Cluster Impact (ECI) 
where a better controlled beam of energetic clusters is deposited on surfaces J57|. Freiburg's group has shown that 
accelerating the clusters leads to improvements in some properties of the films : depositing slow clusters (energy per 
atom 0.1 eV) produces metal films which can be wiped off easily, but accelerating them before deposition (up to 10 eV 
per atom) results in strongly adhering films (58) . MD simulations of cluster deposition [|58| have explained qualitatively 
this behavior : while low energetic clusters tend to pile up on the substrate leaving large cavities, energetic clusters 
lead to a compact film (Fig. ||). It is interesting to note that, even for the highest energies explored in the MD 
simulations (10 eV per atom), no atoms were ejected form the cluster upon impact. The effect of film smoothening 
is only due to the flattening of the cluster when it touches the substrate. Some caution on the interpretation of these 
simulations in needed because of the very short time scales which can be simulated (some ps) . Similar MD simulations 
of the impact of a cluster with a surface at higher energies have also been performed j5^|. Recently, Palmer's group 
has studied the interaction of Ag clusters on graphite for various incident kinetic energies (between 15 and 1500 
eV). They have shown that, for small (Ag^) clusters, the probability of a cluster penetrating the substrate or not 
critically depends on its orientation relative to the substrate. 



2. Low energy clusters 



Another strategy to grow nanostructures with cluster beams consists in depositing low energy particles 
p3| , ^9 48 Gl] |65| . Ideally, by depositing the clusters with low kinetic energies, one would like to conserve the memory 
of the free cluster phase |l3| to form thin films with original properties. Since the kinetic energy is of the order of 10 
eV per cluster [p6fl, i.e. a few meV per atom which is negligible compared to the binding energy of an atom in the 
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cluster, no f ragme ntation of the clusters is expected upon impact on the substrate. Fig. || suggests that the films 
are porous [ p8|]67|| , which is interesting to keep one of the peculiarities of the clusters : their high surface/volume 
ratio which affects all the physical (structural, electronic ^properties as well as the chemical reactivity (catalysis). 
Concerning deposition of carbon clusters, experiments |l3| , po|] as well as simulations [BSj have shown that the carbon 
clusters preserve their identity in the thick film. Another interesting type of nanostructurcd film grown by cluster 
deposition is the growth of cermets by combining a cluster beam with an atomic beam of the encapsulating material 
]69| . The point is that the size of the metallic particles is determined by the incident cluster size and the concentration 
by the ratio of the two fluxes. Then, these two crucial parameters can be varied independently, in contrast to the 
cermets grown from atomic beams and precipitation upon annealing. 

Cluster beams are generated by different techniques : Multiple Expansion Cluster Source (MECS, [pql ), gas- 
aggregation [ ^2| , ^3pl]f70| , |7l[ | . . . All these techniques produce a beam of clusters with a distribution of sizes, with 
a dispersion of about half the mean size. For simplicity, I will always refer to this mean size. In gas-aggregation 
techniques, an atomic vapor obtained from a heated crucible is mixed with an inert gas (usually Ar or He) and the 
two are cooled by adiabatic expansion, resulting in supersaturation and cluster formation. The mean cluster size can 
be monitored by the different source parameters (such as the inert gas pressure) and can be measured by a time of 
flight mass spectrometer. For further experimental details on this technique, see Refs. |13|j72|. To produce clusters 
of refractory materials, a different evaporation technique is needed : laser vaporization [ 13| , |72|]73] |. A plasma created 
by the impact of a laser beam focused on a rod of is thermalized by injection of a high pressure He-pulse (typically, 
3-5 bars during 150 to 300 /is), which permits the cluster growth. The mean cluster size is governed by several 
parameters such as the helium flow, the laser power and the delay time between the laser shot and the helium pulse. 
As a consequence of the pulsed laser shot, the cluster flux reaching the surface is not continuous but chopped. Typical 
values for the chopping parameters are : active portion of the period ~ 100/is and chopping frequency / = 10Hz. 



3. Other approaches 



Alternatively, one can deposit accelerated clusters onto a buffer layer which acts as a "mattress" to dissipate the 
kinetic energy. This layer is then evaporated, which leads to cluster soft-landing onto the substrate (7|j7^]. The 
advantage of this method is that it is possible to select the mass of the ionized clusters before deposition. However, it 
is difficult with this technique to reach high enough deposition rates to grow films in reasonable times. Vitomirov et 
al |76| deposited atoms onto a rare-gas buffer layer : the atoms first clustered on top and within the layer which was 
afterwards evaporated, allowing the clusters to reach the substrate. Finally, deposition of clusters from a Scanning 
Tunneling Microscope (STM) tip has been shown to be possible, both theoretically |77]] and experimentally Q. 



III. MODELS OF PARTICLE DEPOSITION 



I describe in this section simple models which allow to understand the first stages of film growth by low energy 
cluster deposition. These models can also be useful for understanding the growth of islands from atomic beams in the 
submonolayer regime in simple cases, namely (almost) perfect substrates, irreversible aggregation, etc. and they have 
allowed to understand and quantify many aspects of the growth : for a review of analysis of atomic deposition with 
this kind of models, see Refs. [||J39[(7|Q. The models described below are similar to previous models of diffusing 
particles that aggregate, but such "cluster-cluster aggregation" (CCA) models [|lj do not incorporate the possibility 
of continual injection of new particles via deposition, an essential ingredient for thin film growth. 



Given an experimental system (substrate and cluster chemical nature), how can one predict the growth character- 
istics for a given set of parameters (substrate temperature, incoming flux of clusters . . . )? 

A first idea - the "brute-force" approach - would be to run a Molecular Dynamics (MD) simulation with ab-initio 
potentials for the particular system one wants to study. It should be clear that such an approach is bound to fail since 
the calculation time is far too large for present-day computers. Even using empirical potentials (such as Lennard- 
Jones, Embedded Atom or Tight-Binding) will not do because there is an intrinsically large time scale in the growth 
problem : the mean time needed to fill a significant fraction of the substrate with the incident particles. An estimate 
of this time is fixed by tML, the time needed to fill a monolayer : tuL — l/F where F is the particle flux expressed 
in monolayers per second (ML/s). Typically, the experimental values of the flux are lower than lML/s, leading to 
t]\jL > Is. Therefore, there is a time span of about 13 decades between the typical vibration time (10 _13 s, the lower 
time scale for the simulations) and tuL, rendering hopeless any "brute- force" approach. 
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There is a rigorous way |82|] of circumventing this time span problem : the idea is to "coarsen" the description 
by defining elementary processes, an approach somewhat reminiscent of the usual (length, energy) renormalization 
of particle physics [GJ. One "sums up" all the short time processes (typically, atomic thermal vibrations) in effective 
parameters (transition rates) valid for a higher level (longer time) description. I will now briefly describe this rigorous 
approach and then proceed to show how it can be adapted to cluster deposition. 



A. Choosing the elementary processes 

Voter Q showed that the interatomic potential for any system can be translated into a finite set of parameters, 
which then provides the exact dynamic evolution of the system. Recently, the same idea has been applied to Lennard- 
Jones potentials Q by using only two parameters. The point is that these coarse-grained, lattice-gas approach needs 
orders of magnitude less computer power than the MD dynamics described above. One can understand the basic idea 
by the following simple example : for the MD description of the diffusion of an atom by hopping, one has to follow 
in detail its motion at the picosecond scale, where the atom mainly oscillates in the bottom of its potential well. 
Only rarely at this time scale will the atom jump from site to site, which is what one is interested on. Voter showed 
that, provided some conditions are met concerning the separation of these two time scales, and restricting the motion 
to a regular (discrete) lattice (see |8^j for more details), one could replace this "useless" information by an effective 
parameter taking into account all the detailed motion of the atom within the well (including the correlations between 
the motions of the atom and its neighbors) and allowing a rapid evaluation of its diffusion rate. 

Unfortunately, this rigorous approach is not useful for cluster deposition, because the number of atomic degrees of 
freedom (configurations) is too high. Instead, one chooses - from physical intuition - a "reasonable" set of elementary 
processes, whose magnitudes are used as free parameters. This allows to understand the role of each of these elementary 
processes during the growth and then to fit their value from experiments (Fig. |i|). These are the models which I will 



study in this paper, with precise examples of parameter fit (see section VI). Examples of such fits from experimental 
data for atomic deposition include homoepitaxial growth of GaAs(OOl) H], of Pt(100) J85| or several metal(lOO) 
surfaces |8(|. Of course, fitting is not very reliable when there are too many (almost free) parameters. An interesting 
alternative are intermediate cases, where parameters are determined from known potentials but with a simplified 
fitting procedure taking into account what is known experimentally of the system under study : see Ref. |S7| ] for a 
clear example of such a possibility. 



B. Predicting the growth from the selected elementary processes 



To be able to adjust the values of the elementary processes from experiments, one must first predict the growth 
from these processes. The oldest way is to write "rate-equations" which describe in a mean-field way the effect of 
these processes on the number of isolated particles moving on the substrate (called monomers) and islands of a given 
size. The first author to attempt such an approach for growth was Zinsmeister ]88| in 1966, but the general approach 
is similar to the rate-equations first used by Smoluchovsky for particle aggregation |p9[ . In the seventies, many papers 
dealing with better mean-field approximations and applications of these equations to interpret experimental systems 
were published. The reader is referred to the classical reviews by Venables and co-workers [ po|pl| and Stoyanov and 
Kaschiev J9^] for more details on this approach. More recently, there have been two interesting improvements. The 
first is by Villain and co-workers which have simplified enormously the mathematical treatment of the rate-equations, 
allowing one to understand easily the results obtained in a variety of cases [ p3|]94| ]. Pimpinelli et al. have recently 
published a summary of the application of this simplified treatment to many practical situations using a unified 
approach p5[ |. The second improvement is due to Bales and Chrzan p6[ | who have developed a more sophisticated 
self-consistent rate-equations approach which gives better results and allows to justify many of the approximations 
made in the past. However, these analytical approaches are mean-field in nature and cannot reproduce all the 
characteristics of the growth. Two known examples are the island morphology and the island size distribution 96 1 . 

The alternative approach to predict the growth are Kinetic Monte-Carlo (KMC) simulations. KMC simulations 
are an extension of the usual Monte Carlo [p7|-|99]| algorithm and provide a rigorous way of calculating the dynamical 
evolution of a complicated system where a large but finite number of random processes occur at given rates. KMC 
simulations are useful when one chooses to deal with only the slowest degrees of freedom of a system, these variables 
being only weakly coupled to the fast ones, which act as a heat bath |9{j]. The "coarsened" description of film growth 
(basically, diffusion) given above is a good example p^,p6,100 103[], but other applications [Ejil of KMC simulations 



include interdiffusion in alloys, slow phase separations . . . The principle of KMC simulations is straightforward : one 
uses a list of all the possible processes together with their respective rates v pro and generates the time evolution of 
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the system from these processes taking into account the random character of the evolution. For the simple models 
of film growth described below, systems containing up to 4000 x 4000 lattice sites can be simulated in a reasonable 
time (a few hours), which limits the finite size effects usually observed in this kind of simulation. Let me now discuss 
in some detail the way KMC simulations are implemented to reproduce the growth, once a set of processes has been 
defined, with their respective v pro taking arbitrary values or being derived from known potentials. 

There are two main points to discuss here : the physical correctness of the dynamics and the calculation speed. 
Concerning the first point, it should be noted that, originally |)7|], Monte Carlo simulations aimed at the description 
of the equations of state of a system. Then, " the MC method performs a " time" averaging of a model with (often 
artificial) stochastic kinetics [...] : time p lays the role of a label characterizing the sequential order of states, and 
need not be related to the physical times" [ 104 1 . One should be cautious therefore on the precise Monte Carlo scheme 
used for the simu lation wh en attempting at describing t he k i neti cs of a system, as in KMC simulations. For example, 
there are doubts 105 10E] about some simulation work [ 107, 108| carried out using Kawasaki dynamics. This point is 
discussed in great detail in Ref. |106[ . 

Let me address now the important problem of the calculation speed. One could naively think of choosing a time 
interval At smaller than all the relevant times in the problem, and then repeat the following procedure : 

(1) choose one particle randomly 

(2) choose randomly one of the possible processes for this particle 

(3) calculate the probability p pro of this process happening during the time interval At (p pro = v pro At} 

(4) throw a random number p r and compare it with p pro : if p pro < p r perform the process, if not go to the next 
step 

(5) increase the time by At and goto (1) 

This procedure leads to the correct kinetic evolution of the system but might be extremely slow if there is a large 
range of probabilities p pro for the different processes (and therefore some p pro <C 1). The reason is that a significant 
fraction of the loops lea ds to rejected moves, i.e. to no evolution at all of the system. 

Instead, Bortz et al. [10E] have proposed a clever approach to eliminate all the rejected moves and thus reduce 
dramatically the computational times. The point is to choose not the particles but the processes, according to their 
respective rate and the number of possible ways of performing this process (called Q pro ). This procedure can be 
represented schematically as follows : 

(1) update the list of the possible ways of performing the processes VL pro 

(2) randomly select one of the process, weighting the probability of selection by the process rate v pro and Q pro : 

P P ro = {Vprc&pro) I (j2 pr ocesses ^pro^ P ro^j 

(3) randomly select a particle for performing this process 

(4) move the particle 



(5) increase the time by dt = (J2 processes ^ P roV P ro 

(6) goto (1) 



A specific example of such a scheme for cluster deposition is given below (Section III C ). Note that the new 
procedure implies a less intuitive increment of time, and that one has to create (and update) a list of all the £l pro 
constantly, but the acceleration of the calculations is worth the effort. 

A serious limitation of KMC approaches is that one has to assume a finite number of local environments (to 
obtain a finite number of parameters) : this confines KMC approaches to regular lattices, thus preventing a rigorous 
consideration of elastic relaxation, stress effects . . . everything that affects not only the number of first or second nearest 
neighbors but also their precise position. Indeed, considering the precise position as in MD simulations introduces 
a continuous variable and leads to an infinite number of po ssib le configurations or processes. Stress effects can be 
introduced approximately in KMC simulations, for example |11C] by allowing a variation of the bonding energy of an 
atom to an island as a function of the island size (the stress depending on the size), but it is unclear how meaningful 



these approaches a re (s ee also Refs. [ 111 |). I should quote here a recent proposition [ 1 12 1 inspired on the old Frenkel- 



Kontorova model 113 which allows to incorporate some misfit effects in rapid simulations. It remains to explore 
whether such an approach could be adapted to the KMC scheme. 



C. Basic elementary processes for cluster growth 

What is likely to occur when clusters are deposited on a surface ? I will present here the elementary processes which 
will be used in cluster deposition models : deposition, diffusion and evaporation of the clusters and their interaction 
on the surface (Figs. || and ||). The influence of surface defects which could act as traps for the particles is also 
addressed. 



7 



A simple physical rationale for choosing only a limited set of parameters is the following (see Fig. @). For any given 
system, there will be a "hierarchy" of time scales, and the relevant ones for a growth experiment are those much lower 
than £ml — l/F. The others are too slow to act and can be neglected. The hierarchy of time scales (and therefore the 
relevant processes) depend of course on the precise system und er st udy. It should be noted that for cluster deposition 
the situation is somewhat simpler than for atom deposition [114] since many el emen tary processes are very slow. 
For example, diffusion of clusters on top of an already formed island is very low [ 115 1 , cluster detachment from the 
islands is insignificant and edge diffusion is not a n elem entary process at all since the cluster cannot move as an entity 
over the island edge (as I will discuss in section VII B , the equivalent process is cluster-cluster coalescence by atomic 
motion). Let me now discuss in detail each of the elementary processes useful for cluster deposition. 

The first ingredient of the growth, deposition, is quantified by the flux F, i.e. the number of clusters that are 
deposited on the surface per unit area and unit time. The flux is usually uniform in time, but in some experimental 
situations it can be pulsed, i.e. change from a constant value to over a given period. Chopping the flux can affect 
the growth of the film significantly 116], and I will take this into account when needed (Section VI B 3). 

The second ingredient is the diffusion of the clusters which have reached the substrate. I assume that the diffusion 
is brownian, i.e. the particle undergoes a random walk on the substrate. To quantify the diffusion, one can use both 
the usual diffusion coefficient D or the diffusion time r, i.e. the time needed by a cluster to move by one diameter. 
These two quantities are connected by D = d 2 /{Ar) where d is the diameter of the cluster. Experiments show that 
the diffusion coefficient of a cluster can be surprisingly large, comparable to the atomic diffusion coefficients. The 
diffusion is here supposed to occur on a perfect substrate. Real surfaces always present some defects such as steps, 
vacancies or adsorbed chemical impurities. The presence of these defects on the surface could significantly alter the 
diffusion of the particles and therefore the growth of the film. I will include here one simple kind of defect, a perfect 
trap for the clusters which definitively prevents them from moving. 

A third process which could be present in growth is re- evaporation of the clusters from the substrate after a time 
T e . It is useful to define X$ = \/F)T e the mean diffusion length on the substrate before desorption. 

The last simple process I will consider is the interaction between the clusters. The simplest case is when aggregation 
is irreversible and particles simply remain juxtaposed upon contact. This occurs at low temperatures. At higher 
temperatures, cluster-cluster coalescence will be active (Fig. ||). Thermodynamics teaches us that coalescence should 
always happen but without specifying the kinetics. Since many clusters are deposited on the surface per unit time, 
kinetics is here crucial to determine the shape of the islands for med on the substrate. A total comprehension of the 
kinetics is still lacking, for reasons that I will discuss later (Section VII B| ) . I note that the shape of the clusters and the 
islands on the surface need not be perfectly spherical, even in the case of total coalescence. Their interaction with the 
substrate can lead to half spheres or even flatter shapes depending on the contact angle. Contrary to what happens 
for atomic deposition, a cluster touching an island forms a huge number of atom-atom bonds and will not detach from 
it. Thus, models including reversible particle-particle aggregation pl| arc not useful for cluster deposition. 

The specific procedure to perform a rapid KMC simulation of a system (linear size L) when deposition, diffusion 
and evaporation of the monomers are included is the following. The processes are : deposition of a particle (vdepo = F, 
^depo = L 2 (it is possible to deposit a particle on each site of the lattice)), diffusion of a monomer (vdiff — 1/ T i 
^diff = pL 2 where p is the monomer density on the surface) and evaporation of a monomer (v evap — l/r e , Vl evap = 
pL 2 ). For each loop, one calculates two quantities pdrop = F '/ (F ' + + 7)) andp^i/ = (p / T ) / (F + p{^- + -))■ Then, 
one throws a random number p (0 < p < 1) and compare it to Pdrop and pdif- If p < Pdrop, & particle is deposited in a 
random position; if p > Pdrop +Pdif, a monomer (randomly selected) is removed, otherwise we just move a randomly 
chosen monomer. After each of these possibilities, one checks whether an aggregation has taken place (which modifies 
the number of monomers on the surface, and therefore the number of possible diffusion or evaporation moves), increase 
the time by dt — 1/(FL 2 + pL 2 {^- + 7)) and go to the next loop. 



The usual game for theoreticians is to combine these elementary processes and predict the growth of the film. 
However, experimentalists arc interested in the reverse strategy : from (a set of) experimental results, they wish to 
understand which elementary processes are actually present in their growth experiments and what are the magnitudes 
of each of them, what physicists call understanding a phenomenon. The problem, of course, is that with so many 
processes, many combinations will reproduce the same experiments (see specific examples below). Then, some clever 
guesses are needed to first identify which processes are present. For example, if the saturation island density does not 
change when flux (or substrate temperature) is changed, one can guess that nucleation is mostly occurring on defects 
of the surface. 

In view of these difficulties, next section is devoted to predict the growth when the microscopic processes (and 
their values) are known. After, in Section [y|, I propose a detailed procedure to identify and quantify the microscopic 
process from the experiments. Finally, Section v| reviews the experimental results obtained for cluster deposition. 
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IV. PREDICTING GROWTH WITH COMPUTER SIMULATIONS 



The scope of this section is to find formulas or graphs to deduce the values of the microscopic processes (diffusion, 
evaporation . . . ) from the observed experimental quantities (island density, island size histograms . . . ). The "classical" 
studies J90|-p^| have focused on the evolution of the concentration of islands on the surface as a function of time, 
and especially on the saturation island density, i.e. the maximum of the island density observed before reaching a 
continuous film. The reason is of course the double possibility to calculate it from rate-equations and to measure it 
experimentally by conventional microscopy. I will show other interesting quantities such as islan d size dis tributions 



which are measurable experimentally and have been recently calculated by computer simulations [117 122]. 

I will study the two limiting cases of pure juxtaposition and total coalescence (which are similar to two and three 
dimensional growth in atomic deposition terminology) separately. Experimentally, the distinction between the two 
cases can be made by looking at the shape of the supported islands : if they are circular (and larger than the incident 
clusters) they have been formed by total coalescence; if ramified by pure juxtaposition (see several examples below, 



section VI) 



In both cases, I analyze how the growth proceeds when different processes are at work : diffusion, evaporation, 
defects acting as traps, island mobility ... In the simulations, I often take the diffusion time r to be the unit time : in 
this case, the flux is equivalent to the normalized flux <j> (see the Table) and the evaporation time corresponds to r e /r. 
The growth is characterized by the kinetics of island formation, the value of the island concentration at saturation 
N sa t (i.e. the maximum value reached before island-island coalescence becomes important) and the corresponding 
values of the thickness e sa t and condensation coefficient C sa t , useful when evaporation is important (the condensation 
coefficient is the ratio of matter actually present on the substrate over the the total number of particles sent on the 
surface (also called the thickness e = Ft), see Table I). 

I also give the island size distributions corresponding to each growth hy pothesis. The se have proven useful as a tool 
for experimentalists to distinguish between different growth mechanisms [115 121 , 122f| . By size of an island, I mean 



the surface it occupies on the substrate. For "two dimensional" islands (i.e. formed by pure juxtaposition), this is the 
same as the island mass, i.e. its number of monomers. For "three dimensional" islands (formed by total coalescence), 
their projected surface is the easiest quantity to measure by microscopy. It should be noted that for three dimensional 
islands, their projected surface for a given mass depends on their shape, which is assumed here to be pyramidal (close 



to a half-sphere). It has been shown p6,101] that by normalizing the size histograms, one obtains a "universal" size 



distribution independent of the coverage, the flux or the substrate temperature for a large range of their values. 

A. Pure juxtaposition : growth of one cluster thick islands 

I first study the formation of the islands in the limiting case of pure juxtaposition. This is done for several growth 
hypothesis. The rate-equations treatment is given in Appendix A. 

1. Complete condensation 

Let me start with the simplest case where only diffusion takes place on a perfect substrate (no evaporation). Fig 
|§|a shows the evolution of the monomer (i.e. isolated clusters) and island densities as a function of deposition time. 

We see that the monomer density rapidly grows, leading to a rapid increase of island density by monomer- monomer 
encounter on the surface. This goes on until the islands occupy a small fraction of the surface, roughly 0.1% (Fig. 
|jj|a). Then, islands capture efficiently the monomers, whose density decreases. As a consequence, it becomes less 
probable to create more islands, and their number increases more slowly. When the coverage reaches a value close 
to 15% (Fig. ^Jd), coalescence will start to decrease the number of islands. The maximum number of islands at 
saturation N sat is thus reached for coverages around 15%. Concerning the dependence of N sat as a function of the 
model parameters, it has been shown that the maxim um n umber of islands per unit area formed on the surface scales 
as N sat ~ (F/D) 1 / 3 |9p|j92]| . Recent simulations |)6| , |ll5| and theoretical analysis |)3| have shown that the precise 
relation is N sat — 0.53(i 7, r)"- 36 for the ramified islands produced by pure juxtaposition (Fig. |lO| ). 

It should be noted that if cluster diffusion is v anishingly small, the above relation does not hold : instead, film 



such a situation has been given in Ref. [124] 



growth proceeds as in the percolation mode l [ 123 1 , by random paving of the substrate. An experimental example of 
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2. Evaporation 



What happens when evaporation is also included ? Fig ||b shows that now the monomer density becomes roughly a 
constant, since it is now mainly determined by the balancing of deposition and evaporation. As expected, the constant 
concentration equals Fr e (solid line) . Then the number of islands increases linearly with time (the island creation rate 
is roughly proportional to the square monomer concentration, see Appendix A). One can also notice that only a small 
fraction (1/100) of the monomers do effectively remain on the substrate, as shown by the low condensation coefficient 
value at early times. This can be understood by noting that the islands grow by capturing only the monomers 
that are deposited within their "capture zone" (comprised between two circles of radius R and R + Xs)- The other 
monomers evaporate before reaching the islands. When the islands occupy a significant fraction of the surface, they 
capture rapidly the monomers. This has two effects : the monomer density starts to decrease, and the condensation 
coefficient starts to increase. Shortly after, the island density saturates and starts to decrease because of island-island 
coalescence. Fig. ^ shows the evolution of the maximum island density in the presence of evaporation. A detailed 
analysis of the effect of monomer evaporation on the growth is given in Ref. [117|, where is also discussed the regime 
of " direct impingement" which arises when Xs < 1 : islands are formed by direct impingement of incident clusters as 
first neighbors of previously adsorbed clusters, and grow by direct impingement of clusters on the island boundary. 
A summary of the results obtained in the various regimes spanned as the evaporation time r e decreases is given in 
Appendix A. 



3. Defects 



I briefly treat now the influence of a very simple kind of defect : a perfect trap for the diffusing particles. If a 
particle enters such a defect site, it becomes trapped at this site forever. If such defects are present on the surface 
[125 1 they will affect the growth of the film only if their number is higher than the number of islands that would have 
been created without defects (for the same values of the parameters). If this is indeed the case, monomers will be 
trapped by the defects at the very beginning of the growth and the number of islands equates the number of defects, 
whatever the diffusivity of the particles. The kinetics of island formation is dramatically affected by the presence of 
defects, the saturation density being reached almost immediately (Fig. [Tl|) . 



4- Island mobility 



The consequences of the mobility of small islands have not received much attention. One reason is that it is difficult 
to include island mobility in rate equations treatments. A different (though related!) reason is that (atomic) islands 
are expected to be almost immobile in most homoepitaxial systems. However, several studies have shown the following 
consequences of island mobility for the pure juxtaposition case and in the absence of evaporation. First, the saturation 
island density is changed p^ JlO^ , p^ p8| : one obtains N sat = 0.3(F/D) A2 (Fig. [[(]) if all islands are mobile, with 
a mobility inversely proportional to their size [ 102 1 . Second, the saturation island density is reached for very low 
coverages (Fig. [ft] and Ref. JTolfl). This can be explained by a dynamical equilibrium between island formation and 
coalescence taking place at low coverages thanks to island diffusion. If only monomers a re ab le to move, islands can 
coalesce (static coalescence) only when the coverage is high enough (roughly 10-15%, (92,117]). Then, the saturation 
island density is reached in this case for those coverages. Instead, when islands can move, the so calle d dy namical 
coalescence starts from the beginning of the growth and the balance is e stablished at very low coverages [ 102 1 . Third, 
the island size distribution is sharpened by the mobility of the islands [ 118, 120, 127 , 128 1 . To my knowledge, there is 
no prediction concerning the growth of films with evaporation when islands are mobile. 



5. Island size distributions 



Fig [l^ shows the evolution of the r eseated 101 121 ] island size distributions as a function of the evaporation time for 
islands formed by juxtaposition [117]. Size distributions are normalized by the mean island size in the following way : 
one defines p(s/s m ) = n s /Nt as the probability that a randomly chosen island has a surface s when the average surface 
per island is s m — 9/N t , where n s stands for the number of islands of surface s, N t is the total number of islands and 
8 for the coverage of the surface. It is clear that the distributions are significantly affected by the evaporation, smaller 
islands becoming more numerous when evaporation increases. This trend can be qualitatively understood by noting 
that new islands are created continuously when evaporation is present, while nucleation rapidly becomes negligible 
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in the complete condensation regime. The reason is that islands are created (spatially) homogeneously in the last 
case, because the positions of the islands are correlated (through monomer diffusion), leaving virtually no room for 
further nucleation once a small portion of the surface is covered (9 ~ 0.05). In the limit of strong evaporation, islands 
are nucleated randomly on the surface, the fluctuations leaving large regions of the surface uncovered. These large 
regions can host new islands even for relatively large coverages, which explains that there is a large proportion of 
small (s < s m ) islands in this regime. 



B. Total coalescence : growth of three dimensional islands 

If clusters coalesce when touching, the results are slightly different from those given in the preceding section, mainly 
because the islands occupy a smaller portion of the substrate at a given thickness. Therefore, in the case of complete 
condensation for example, saturation arises at a higher thickness (Fig. |ll| ) even if the coverage is approximately 
the same (matter is "wasted" in the dimension perpendicular to the substrate). However, the main qualitative 
characteristics of the growth correspond to those detailed in the preceding section. Fig. |l^ shows the evolution of 
the maximum island density in that case, where the three-dimensional islands are assumed to be roughly half spheres 
(actually, pyramids were used in these simulations which were originally intended for atomic deposition fl29f |). The 
analytical results obtained from a rate-equations treatment are given in Appendix B. If the islands are more spherical 
(i.e. the contact angle is higher), a simple way to adapt these results on the kinetic evolution of island concentration 
(Fig. [ll]) is to multiply the thickness by the appropriate form factor, 2 for a sphere for example. Indeed, if islands 
are spherical, the same coverage is obtained for a thickness double than that obtained for the case of half-spheres 
(there are two identical half spheres). This is a slight approximation since one has to assume that the capture cross 



section (which governs the growth) is identical for the two shapes : this is not exactly true [117] but is a very good 
approximation. 

Fig [0| shows the evolution of the rescaled island size distributions for three dimensional islands (pyramids) in 
presence of evaporation. I recall that size means here the projected surface of the island, a quantity which can be 
measured easily by electron microscopy. We note the same trends as for the pure juxtaposition case. 

Fig |l5] shows the evolution of the rescaled island size distributions for pyramidal islands nucleating on defects. Two 
main differences can be noted. First, the histograms are significantly narrower than in the preceding case, as had 



already been noted in experimental studies [ 130 1 . This can be understood by noting that all islands are nucleated at 
almost the same time (at the very beginning of growth). The second point is that the size distributions are sensitive 
to the actual coverage of the substrate, in contrast with previous cases. In other words, there is no perfect rescaling 
of the data obtained at different coverages, even if rescaling for different fluxes or diffusion times has been checked. 



C. Other growth situations 

I briefly address in this paragraph other processes which have not been analyzed here. A possible (but difficult 
to study) process is a long range interaction between particles (electrostatical or throu gh the substrate). There is 
some experimental evidence of this kind of interaction for the system Au/KCl(100) ]131J but to my knowledge, it has 
never been incorporated in growth models. Chemical impurities adsorbed on the substr ate c an change the growth in 



conventional vacuum, and these effects are extremely difficult to understand and control [132]. Of course, many other 
possible processes have not been addressed in this review, such as the influence of strain, of extended defects as steps 
or vacancy islands . . . 



V. HOW TO ANALYZE EXPERIMENTAL DATA 

Figures [T^, [l3| constitute in some sense "abacuses" from which one can determine the value of the microscopic 
parameters (diffusion, evaporation) if the saturation island density is known. The problem is : does the measured 
island density correspond to the defect concentration of the surface or to homogeneous nucleation? If the latter is 
true : which curve should be used to interpret the data? In other words, is evaporation present in the experiments 
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and what is the magnitude of r e ? I will now give some tricks to first find out which processes are relevant and then 
how they can be quantified. 

Let's concentrate first on the presence of defects. One possibility is to look at the evolution of N sat with the flux. 
As already explained, if this leaves unaffected the saturation density, nucleation is occurring on defects. A similar 
test can be performed by changing the substrate temperature, but there is the nagging possibility that this changes 
the defect concentration on the surface. It is also possible to study the kinetics of island nucleation, i.e. look at the 
island concentration as a function of thickness or coverage. The presence of defects can be detected by the fact that 
the maximum island density is reached at very low coverages (typically less than 1%, see Fig. pd| ) and/or by the 
fact that the nucleation rate (i.e. the derivative of the island density) scales as the flux and not as the square flux 
: see Section 3 of Ref. |)2| for more details. One should be careful however to check that all the islands, even t hose 
containing a few particles, are visible in the microscope images. This is a delicate point for atomic deposition [133] 
but should be less restrictive for clusters since each cluster has already a diameter typically larger than a nanometer. 
Of course all this discussion assumes that the defects are of the "ideal" kind studied here, i.e. perfect traps. If atoms 
can escape from the defects after some time, the situation is changed but I am unaware of studies on this question. 

The question of evaporation is more delicate. First, one should check whether particle reevaporation is important. 
In principle, this can be done by measuring the condensation coefficient, i.e. the amount of matter present on the 
surface as a function of the amount of matter brought by the beam. If possible, this measure leaves no ambiguity. 
Otherwise, the kinetics of island creation is helpful. If the saturation is reached at low thicknesses (e sa t < -5 ML), 
this means that evaporation is not important. Another way of detecting particle evaporation is by studying the 
evolution of the saturation island density with the flux : in the case of 2D growth (Fig. [l(]), the expo nent is 0.36 
when evaporation is negligible but roughly 0.66 when evaporation significantly affects the growth |117|. There are 
similar differences for 3d islands : the exponent changes from 0.29 to 0.66 |134| (Fig. Il3f). Suppose now that one 
finds that evaporation is indeed important : before being able to use Fig. nfl or Fig. H3, one has to know the 
precise value of r e . One way to find out is to make a precise fit of the kinetic evolution of the island density or the 
condensation coefficient (see Section VI B 2 for an example). In next paragraph, I show how to find r e if one knows 
only the saturation values of the island density and the thickness. 



As a summary, here is a possible experimental strategy to analyze the growth. First, get a series of micrographs 
of submonolaycr films as a function of the thickness. The distinction between the pure juxtaposition and total 
coalescence cases can be easily made by comparing the size of the supported islands to the (supposedly known) size 
of the incident clusters. Also, if the islands are spherical, this means that coalescence has taken place, if they are 
ramified that clusters only juxtapose upon contact. Of course, all the intermediate cases are possible (see the case of 
gold clusters below). One can calculate the ratio of deposited thickness over the coverage : if this ratio is close to 1, 
islands are flat (i.e. one cluster thick), otherwise three dimensional (unless there is evaporation). 

From these micrographs, it is possible to measure the island density as a function of the thickness. Fig. |ll| 
should be now helpful to distinguish between the different growth mechanisms. For example, if the saturation island 
density is obtained for large thicknesses (typically more than 1ML), then evaporation is certainly relevant and trying 
to measure the condensation coefficient is important to confirm this point. It is clear from Figs. [L(j H that the 
knowledge of N sat alone cannot determine r e since many values of r and r e can lead to the same N sat . In the 2D 
case, the values of the microscopic parameters can be obtained by noting that the higher the evaporation rate, the 
higher the amount of matter "wasted" for film growth (i.e. re-evaporated). One therefore expects that the smaller 
r e , the higher e sa t, which is confirmed by Fig. |l6|a. Therefore, from the (known) value of e sa t, one can determine the 
value of the evaporation parameter r\ = FrXg (Fig- |l6|a). Once rj is known, X$ is determined from Fig. [l6[ b since 
N sat is known. Ft can afterwards be determined (from X$ and rf). This is only valid for X$ ^S> 1 117 ], a condition 
always fulfilled in experiments. 

The 3d case is more difficult since the same strategy (measuring N sat and e sa t) fails. The reason is that in the 

1/2 

limit of high evaporation, e sat goes as e sa t ~ N sat 1 , thus bringing no independent information on the parameters 
[ 134 1 . The same is true for the condensation coefficient at saturation C sa t, which is a constant, i.e. independent of 



the value of r e or the normalized flux (see Fig. 17b). This counterintuitive result (one would think that the higher 



the evaporation rate, the smaller the condensation coefficient at saturatio n) c an be understood by noting that in this 
limit, islands only grow by direct impingement of particles within them [134] and therefore X$ (or r e ) has no effect 
on the growth. Fortunately, in many experimental situations, the limit of high evaporation is not reached and one 
"benefits" from (mathematical) crossover regimes where these quantities do depend on the precise values of r e . Figs. 
|l7| give the evolutions of C sa t and e sa t as a function of N sat for different values of r e and F. Then, knowing e sa t and 
N sat leads to an estimation of r e from Fig. |l7| a which can be confirmed with Fig. |l7j b provided C sa t is known. 
To conclude, let me note that a saturation thickness much smaller than 1ML can also be attributed to island 
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mobility. This is a subtle process and it is difficult to obtain any information on its importance. We note that 
interpreting data as not affected by island diffusion when it is actually present leads to errors on diffusion coefficients 
of one order of magnitude or more depending on the value of Ft (see Fig. [lO]). Finally, one should be careful 
in interpreting the N t vs. thickness curves since most observations are not made in real time (as in the computer 



simulations) and there can be post-deposition evolutions (see for example Ref. [ 135 1 for such complications in the case 
of atomic deposition) . 



VI. EXPERIMENTAL RESULTS 



I review in this section the experimental results obtained these last years for low-energy cluster deposition, mainly 
in the submonolayer regime. The scope is double : first, 1 want to give some examples on how to analyze experiments 
(as indicated in Section |V|) and second, I will show that from a comparison of experiments and models one can deduce 
important physical quantities characterizing the interaction of a cluster with a surface (cluster diffusivity) and with 
another cluster (coalescence). The following can be read with profit by those interested only in atomic deposition 
as examples of interpretation since these elementary processes are relevant for some cases of atomic deposition. One 
should be careful that some mechanisms which are specific to atomic deposition (transient mobility, funnelling, . . . ) 
are not discussed here (see [136]). Also, g row th without cluster diffusion has to be interpreted in the framework of 
the percolation model as indicated above [124]. 

Before analyzing experimental data, it is important to know how to make the connection between the units used 
in the programs and the experimental ones (see also Table I). In the program, the unit length is the diameter of a 
cluster. In the experiments, it is therefore convenient to use as a surface unit the site, which is the projected surface 
of a cluster 7re? 2 /4 where d is the mean incident cluster diameter. The flux is then expressed as the number of clusters 
reaching the surface per second per site (which is the same as ML/s) and the island density is given per site. The 
thickness is usually computed in cluster monolayers (ML), obtained by multiplying the flux by the deposition time. 
The coverage - the ratio of the area covered by the supported islands over the total area - has to be measured on the 
micrographs. 



A. A simple case : 562300 clusters on graphite HOPG 



I start with the case of antimony clusters c onta ining 2300 (±600) atoms deposited on graphite HOPG since here 
the growth has been thoroughly investigated ]115| |. I first briefly present the experimental procedure and then the 
results and their interpretation in terms of elementary processes. 



1. Experimental procedure 

As suggested in the preceding Section, various samples are prepared for several film thicknesses, incident fluxes and 
the substrate temperatures. For films grown on Highly Oriented Pyrolitic Graphite (HOPG), before deposition at 
room temperature, freshly cleaved graphite samples are annealed at 500°C during 5 hours in the deposition chamber 
(where the pressure is ~ 10~ 7 Torr) in order to clean the surface. The main advantage of HOPG graphite conveniently 
annealed is that its surfaces consist mainly of defect-free large terraces (~ 1/im) between steps. It is also relatively 
easy to observe these surfaces by electron or tunneling microscopy |115| ]. Therefore, deposition on graphite HOPG 
is a good choice to illustrate the interplay between the different elementary processes which combine to lead to the 
growth. After transfer in air, the films are observed by Transmission Electron Microscopy (TEM) (with JEOL 200CX 
or TOP CON electron microscopes operating at 100 kV in order to improve the contrast of the micrographs). 



2. Results 



Fig. 18 1 shows a general vie w of the morphology of the antimony submonolayer film for e = 0.14 ML and 
T s = 353K. A detailed analysis [115] of this kind of micrographs shows that the ramified islands are formed by 
the juxtaposition of particles which have the same size distribution as the free clusters of the beam. From this, we 
can infer two important results. First, clusters do not fragment upon landing on the substrate as indicated in the 
introduction. Second, antimony clusters remain juxtaposed upon contact and do not coalesce to form larger particles 
(option (a) of Fig. ^). 
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From a qualitative point of view, Fig. |Ts| a also shows that the clusters are able to move on the surface. Indeed, 
since the free clusters are deposited at random positions on the substrate, it is clear that, in order to explain the 
aggregation of the clusters in those ramified islands, one has to admit that the clusters move on the surface. How 
can this motion be quantified? Can we admit that diffusion and pure juxtaposition are the only important physical 
phenomena at work here? 

Fig. |Tg| a shows the evolution of the island density as a function of the deposited thickness. We see that the 
saturation island density N sa t is reached for e ~ 0.15 AIL. This indicates that evaporation or island diffusion are 
not important in this case. Therefore, we guess that the growth should be described by a simple combination of 
deposition, diffusion of the incident clusters and juxtaposition. Th is has been confirmed in several ways. I only 
give three different confirmations, directing the reader to Ref. |l!5|| for further details. First, a comparison of the 
experimental morphology and that predicted by models including only deposition, diffusion and pure juxtaposition 
shows a very good agreement (Fig. [i"8|b). Second, Fig. |l9| b shows that the saturation island density accurately 
follows the prediction of the model when the flux is varied. I recall that if the islands were nucleated on defects of 
the surface, the density would not be significantly affected by the flux. 

Having carefully checked that the experiments are well described by the simple DDA model, I can confidently use 



Fig. |10| to quantify the diffusion of the clusters. As detailed in Ref. [115], one first measures the saturation island 
density for different substrate temperatures. The normalized fluxes (Ft) are obtained from Fig. [lQ[ Knowing the 
experimental fluxes, one can derive the diffusion times and coefficients. The result is a surprisingly high mobility of 
*S'62300 on graphite, with diffusion coefficients of the same order of magnitude as the atomic ones, i.e. 10~ 8 cm 2 s~ 1 
(Fig. 0c). 

The magnitude of the diffusion coefficient is so high that we wondered whether there was any problem in the 
interpretation of the data, in despite of the very good agreement between experiments and growth models described 
above. For example, one could think of a linear diffusion of the incoming clusters, induced by the incident kinetic 
energy of the cluster in the beam (the cluster could "slide" on the graphite surface). This seems unrealistic for two 
reasons : first, in order to explain the low island density obtained in the experiments (see above), it should be assumed 
that the cluster, which has a low kinetic energy (less than 10 eV), can travel at least several thousand nanometers 
before being stopped by friction with the substrate. This would imply that the diffusion is just barely influenced by 
the substrate, which only slows down the cluster. In this case, it is difficult to explain the large changes observed in 
the island density when the substrate temperature varies. Second, we have deposited antimony clusters on a graphite 
substrate tilted to 30° from its usual position (i.e. perpendicular to the beam axis). Then, a linear diffusion of 
the antimony clusters arising from their incident kinetic energy would lead to anisotropic islands (they would grow 



differently in the direction of tilt and its perpendicular). Experiments [HE] show that there is no difference between 
usual and tilted deposits. Therefore we can confidently believe that S'&2300 clusters perform a very rapid Brownian 
motion on graphite surfaces. A similar study has been carried out for S^so on graphite, showing the same order of 



be presented in section VII 



magnitude for the mob ility of the clusters 115]. The microscopic mechanisms that could explain such a motion will 



B. Other experiments 



In this subsection, I try to analyze data obtained in previous studies [32 137 . I provide possible (i.e. not in 
contradiction with any of the data) explanations, with the respective values of the microscopic processes. I stress that 
the scope here is not to make precise fits of the data, but rather to identify the microscopic processes at work and 
obtain good guesses about their respective values. 

1. Slightly accelerated Agieo clusters on HOPG 

Palmer's group |32| has investigated the growth of films by Agi^o cluster deposition. Fig. |2(] shows the ramified 
morphology of a submonolayer deposit. Although no precise fit is possible given the limited experimental data, the 
island density and size shows that Agi§o clusters are mobile on HOPG. 

2. Sb 36 on a-C 

Small antimony clusters are able to move on amorphous carbon, as demonstrated by Figs, ^l], and by the fact that 
the films are dramatically affected by changing the incident flux |137[ . 
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Fig. pl| a shows that these small clusters gather in large islands and coalesce upon contact. The island density is 
shown in Fig. |2l|b. The maximum is reached for a very high thickness (e ~ 1.8ML), which can only be explained 
by supposing that there is significant reevaporation of 5*636 clusters from the surface. Evaporation of small antimony 
clusters (Sb„ with n < 100) from a-C substrates has also been suggested by other authors Jl^]. A fit, using r e = 20 
deduced from Fig. [l7| a gives, with Ft — 10 -5 for F = 6 10~ 3 clusters szte _1 s _1 , leading to r ~ 2 lO~ 3 s,D = 
2 10~ 12 cm 2 /s, r e = .04s and X$ ~ 6 nm before evaporation, and a condensation coefficient of .2 when the maximum 
island density is reached. However, some authors have argued jl8| that the condensation coefficient is not so low. 
It is interesting to try a different fit of the data - in better agreement with this indication - to give an idea of the 
uncertainties of the fits. For this, I assume that the deposited islands are spherical (solid line of Fig. pTjb) by the 
procedure described in Section [TV]. Here I have taken Ft = 3 10~ 6 /or/ = 6 10~ 3 clusters site~ 1 s~ 1 : leading to 
r ~ 5 10 a, r e = .04s corresponding to D = 8 10~ 12 cm 2 /s, and Xg ~ 11 nm before evaporation, and a condensation 
coefficient of .5 when the maximum island density is reached. Note that the condensation coefficient is, as expected, 
higher than in the previous fit and that the agreement with the experimental island densities for the lowest thicknesses 
is better. Comparing the two fits, it can be seen that the difference in the diffusion coefficient is a factor of 4, and 
a factor 2 in the Xs . This means that the orders of magnitude of the values for the microscopic mechanisms can be 
trusted despite lac k of comprehensive experimental investigation. 

Similar studies [188] have allowed to obtain the diffusion and evaporation characteristic times for other clusters 
deposited on amorphous carbon. For Bigo, one finds D ~ 3 10 cm s and Xs ~ 8 nm and for Inioo ■ D ~ 
4 10 _15 cm 2 s _1 and strong coalescence (the incident clusters are liquid). 



3. A11250 on graphite 

Fig. shows the morphology of a gold submonolayer film obtained by deposition of A11250 (±100) clusters prepared 
by a laser source on graphite in a UHV chamber for different substrate temperatures. 

The structures are strikingly similar to those obtained in the 5&2300 case : large, ramified islands. We can conclude 
that AU250 clusters do move on graphite, and that they do not completely coalesce. A more careful examination of 
the island morphology indicates that the size of the branches is not the same as the size of the incident clusters, as 
was the case for 5*62300- Here the branches are larger, meaning that there is a partial coalescence, limited by the 
kinetics of the growth. This is a very interesting experimental test for coalescence models that are presented later. 
I first try to estimate the diffusion coefficient of the gold clusters. We have to be careful here because the incident 
flux is chopped with the laser frequency, roughly 10Hz. The active portion of the period (i.e. when the flux is "on") 
is ~ 100/xs. 

An analysis of the growth in presence of a chopped flux has been reported elsewhere 
values of N sat as a function of the diffusion time t in these experimental conditions (Fi—6 ML/s) for two hypothesis 
: only the monomers move or islands up to the pentamer move too (island mobility is supposed to be inversely 
proportional to its mass). Note that there is a range of diffusion times (up to two orders of magn itude) which lead to 
the same island saturation value, a strange situation in homogeneous nucleation : see Refs. prq , 

Given the experimental island densities, the diffusion coefficients in both hypothesis are shown in Fig. |23Jd. The 
values of the diffusion coefficient seem too high, especially in the case of exclusive monomer diffusion, but there 
is no experimental evidence of island mobility for the moment. I note however that since the incident clusters do 
significantly coalesce, it is not unreasonable to assume that the smallest islands (which are spherical as the incident 
clusters) can move too. We are presently carrying additional tests (on cluster reevaporation or non brownian cluster 
diffusion) to confirm the observation of such high diffusion coefficients. 



,[138 1 . Fig. ||a shows the 



13§ for details. 



4. A«25o on NaCl 



Given the surprising high mobility of A11250 (±100) on HOPG, it was worth testing gold cluster mobility on other 
substrates. I present here recent results obtained by depositing A1L250 clusters on NaCl [13£]. The high island density 
(Fig. p4|a) shows that gold clusters are not very mobile on this substrate, with an upper limit on t he d iffusion 
coefficient of D ~ 10 _15 cm 2 /s. This is in agreement with the low mobilities observed by other authors [1401 in the 
seventies. The diffraction pattern (Fig. p4|b) is similar to that obtained in Figs. 15 c and d of Ref. |140| . The 
authors interpreted their results with the presence of multi-twinned Au particles with two epitaxial orientations : 
Au(lll)/NaCl(100) and Au(100)/NaCl(100). This is reasonable taking into account the interatomic distances for 



these orientations : d 



Auau(III) 



=0.289nm, d 



Na c l(100) : 



=0.282nm and d 



Au A u(lQ0) : 



=0.408nm l/2d 



JVa c i(100) : 



=0.398nm 



(along the face diagonal). These preliminary results suggest that epitaxy may prevent clusters from moving rapidly 
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on a surface, a result which has also be observed by other groups (see next section). They also show that, at least in 
this case, forming the clusters on the surface by atomic aggregation or depositing preformed clusters does not change 
the orientation nor the diffusion of the clusters on the surface. Work is in progress to determine the prec ise atomic 
structure of the clusters, their orientation on the substrate and their diffusion at higher temperatures 139 1. 



VII. TOWARDS A PICTURE OF CLUSTER DIFFUSION AND COALESCENCE AT THE ATOMIC 

SCALE 



In the preceding sections I have tried to analyze the growth with the help of two main ingredients : diffusion 
of the clusters on the surface and their interaction. I have taken the diffusion as just one number quantifying the 
cluster motion, without worrying about the microscopic mechanisms which could explain it. For atomic diffusion, 
these mechanisms have been extensively studied |28 
compact (111) flat surfaces, diffusion occurs by site to site jumps over bridge sites (the transition state). Therefore, 
diffusion is an activated process and plotting the diffusion constant vs. the temperature yields the height of the 
barrier, which gives information about the microscopies of diffusion. This kind of simple interpretation is not valid 
for cluster diffusion. It is always possible to infer an " activation" energy from an Arrhenius plot (see Fig. [Ts|c) but 
the meaning of this energy is not clear since the precise microscopic diffusion mechanism is unknown. 

Similarly, cluster-cluster coalescence (Fig. ||) has been supposed to be total or null (i.e. pure juxtaposition) but 
without considering the kinetics nor the intermediate cases which can arise (see the experimental results for gold on 
graphite for example). 

In this section, I describe some preliminary results which can shed some light on the microscopic mechanisms leading 
to cluster diffusion or coalescence. 



,136, 141 1 and are relatively well-know. In the (simplest) case of 



A. Diffusion of the clusters 



Before turning to the possible microscopic mechanisms, one must investigate whether cluster diffusion is indeed 
such a general phenomenon. Let me review now the available experimental data concerning the diffusion of three- 
dimensional (3D) clusters. I already presented in the previous section several examples of high cluster mobilities 
over surfaces. In the case of 6*62300 on graphite, mobilities as high as D = 10~ 8 cto 2 s _1 are obtained at room 
temperature, and similar values can be inferred for Ag cluster deposition [ ]6l| , p2| . On a-C substrates, diffusion is not 
that ra pid, but has to be taken into account to understand the growth. More than twenty years ago, the Marseille 
group [ [131 140 142.143] carefully studied the mobility of nanometer-size gold crystallites on ionic substrates (MgO, 
By three different methods, they proved that these 3D clusters - grown by atomic deposition at room 
- are significantly mobile at moderately high te mpe ratures (T ~ 350K). The three different metho ds w ere 



KC1, NaCl). 
temperature 



direct observation under the el ectr on microscope beam [142 

[pl| be i 
lod |13lj| 



the radial distribution functions 
I will focus here on the last method 



comparison of abrupt concentration profiles [14C[ or 
before and after annealing. All these results are carefully reviewed in Ref. [ 143 1 . 



Fig. 



I shows the radial distribution functions of the gold clusters obtained 
just after deposition (the flat curve) and after annealing (the oscillating curve) a similar deposit for a few minutes at 
350K (Fig. 4 of Ref. |l3lj ) . The flat curve is a standard as-grown radial distribution function (see for example Ref. 
Jl0l| ). The other curve is significantly different from the first, although the cluster size distribution remains identical 
(Fig. ^5|). This shows that gold clusters move as an entity on KCl(lOO) at 350K, since the conservation of the size 
distribution rules out atomic exchange between islands (the (EC) mechanism presented below (Section VIIAl). From 
the shape of the radial distribution function some features of the cluster-cluster inte raction c ould be derived, m ainly 
that it is a repulsive interaction. The detailed interaction mechanisms are not clear [ 131 , 143 1 . A different study [ 14C ] 
showed that the clusters were mobile only for a limited amount of time (several minutes), and then stopped. It turns 
out that clusters stop as soon as they reach epitaxial orientation on the substrate. Indeed, the gold(lll) planes can 
orient on the KCl(lOO) surface, reaching a stabl e, minim um energy configuration (for more details on the epitaxial 
orientations of gold clusters on NaCl, see Refs. [ 144, 145|). Therefore, 3D cluster diffusion might be quite a common 
phenomenon, at least when there is no epitaxy between the clusters and the substrate. 

What are the possible microscopic mechanisms? Unfortunately for the field of cluster deposition, recent theoretical 
and experimental work has focused mainly on one atom thick, two-dimensional islands, whose diffusion mechanisms 
might be different from those of 3D islands. The focus on 2D islands is due to the technological impetus provided by 
applications of atomic deposition - notably MBE for which one wants to achieve flat layers. Let's briefly review the 
current state of the understanding of 2D island diffusion to see what inspiration we can draw for 3D cluster diffusion. 
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1. 2D island diffusion mechanisms 



There are two main types of mechanisms proposed to account for 2D island diffusion : single adatom motion and 
collective (simultaneous) atom motion. It should be noted that small islands (less than ~ 15 atoms) ar e likely to 
move by specific mechanisms, depending on the details of the island geometry and atomic energy barriers [ 146 - 148 1 . 
Therefore I concentrate here on larger 2D islands. 



a. Individual mechanisms The most common mechanism invoked to account for 2D island diffusion has been that 
of individual atomic motion. By individual I mean that the movement of the whole island can be decomposed in the 
motion of uncorrelated single atom moves. There are two main examples of such a diffusion : evaporation-condensation 
(EC) and periphery-diffusion (PD). Theoretical investigations on these individual mechanisms have generated much 
interest since it was conjectured that this diffusion constant Di nc i is proportional to the island number of atoms 
(island mass) to some power which depends on the precise mechanism (EC or PD) causing island diffusion but not on 
temperature or the chemical nature of the system. If true, this conjecture would prove very useful, for it would allow 
to determine experimentally the mechanism causing island migration by measuring the exponent and some details of 
the atom diffusion energetics by measuring how Di n d depends on temperature. Unfortunately, recent studies have 
shown that this prediction is too simplistic, as I show now for the two different mechanisms. 

(i) Periphery diffusion 

Fig. Eq shows the el eme ntary mechanism leading to island diffusion via atomic motion on the edge of the island 
(label PD). Assuming [|o|] that 

- each atomic jump displaces the center of mass of the island by a distance of order 1 /N (where N is the number of 
atoms of the island), 

- each edge atom (density n s ) jumps with a rate k ~ exp(—E ec i/kBT) where E ec i is the activation energy for jumping 
from site to site alo ng the border and fc^ is Boltzmann constant 

One obtains [ 105 1 : 



D md ~ kn s l/N 2 ~ exp{-E ed /k B T)N- z ' 2 



(7.1) 



if one postulates that n s , the mean concentration of edge atoms is proportional to the perimeter of the island (i.e. 
to N 1 / 2 ). This equation allows in principle to determine the edge activation energy by measuring the temperature 
dependence of 1)^. 



7.1 



However, recent experiments | 15C | and Kinetic Monte-Carlo simulations |p2,105,149,151| have suggested that Eq 



is wrong. First, the size exponent is not universal but depends on the precise energy barriers for atomic motion (and 
therefore on the chemical nature of the material) and, second, the measured activation energy does not correspond 
to the atomic edge diffusion energy. The point is that the limiting mechanism for island dif fusio n is corner breaking, 
for islands would not move over long distances simply by edge diffusion of the outer atoms [105|. Further studies are 
needed to fully understand and quantify the PD mechanism. 

(ii) Evaporation-Condensation 



An alternative route to diffusion is by exch ange of atoms between the island and a 2D atomic gas. This is the usual 

Atoms can randomly evaporate from the island and atoms belonging to 



152] 



mechanism leading to Ostwald ripening 
the 2D gas can condensate on it (Fig. E( 



This leads to fluctuations in the position of the island center of mass, which 
are difficult to quantify because of the possible correlations in the atomic evaporation and condensation. Indeed, an 
atom which has just evaporated form an island is lik ely to condensate on it again, which cannot be accounted by 
a mean- field theory of island-gas exchange of atoms [ 153 1 . The latter leads to a diffusion coefficient scaling as the 
inverse radius of the island [154 ], while correlations cause a slowing down of diffusion, which scales as the inverse 
square radius of the island [108,151 153 . 



Experimentally, Wen et al. [ |154[ have observed by STM the movement of Ag 2D islands on Ag(100) surfaces. They 
measured a diffusivity almost independent of the island size, which rules out the PD mechanism and roughly agrees 
with their [ 154 1 calculation of the size dependence of the EC mechanism. Since this calculation has been shown to be 



only approximate, further theoretic al an d experimental work is needed to clarify the role of EC in 2D island diffusion. 
However, the work by Wen et al. |l54| has convincingly shown that these islands move significantly and that, for 
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silver, island diffusion is the main route to the evolution of the island size distribution, contrary to what was usually 
assumed (Ostwald ripening exclusively due to atom exchange between islands, via atom diffusion on the substrate). 



b. Collective diffusion mechanisms These indivi dual mechanisms lead in general to relatively slow diffusion of the 
islands (of order 10 _17 cm 2 /s at room temperature [154])- For small clusters, different (and fa ster) mechanisms such 
as dimer sheari ng, i nvolving the simultaneous displacement of a dimer, have been proposed [ |155| . More generally, 
Hamilton et al. [156] have proposed a different mechanism, also involving collective motions of the atoms, which leads 
to fast island motion. By collective I mean that island motion is due to a simultaneous (correlated) motion of (at 
least) several atoms of the islan d. 

Specifically, Hamilton et al. |156{ ] proposed that dislocation motion could cause rapid diffusion of relatively small 
(5 to 50 atoms) homoepitaxial islands on fcc(lll) surfaces. Fig. |7| shows the basic idea : a row of atoms move 
simultaneously from fee to hep sites, thus allowing the motion of the dislocation and consequently of the island center 
of mass. Alternative possibilities suggested by Hamilton et al. for dislocation mediated island motion are the "kink" 
mechanism (the same atomic row moves by sequential but correlated atomic motion) or the "gliding" mechanism 
studied below, where all the a toms of the island move simultaneously. Molecular Dynamics simulations, together with 
a simple analytical approach []l 56f| suggest that for the smallest islands (N < 20) the gliding mechanism is favored, 
for intermediate sizes (20 < N < 100) the dislocation motion has the lowest activation energy, while for the largest 
studied islands (N > 100) the preferential mechanism is that of "kink" dislocati on m otion. It is interesting to quote at 
this point recent direct observations of cluster motion by field ion microscopy [157]. Fig. p8| shows successive images 
of a compact Ir\g cluster moving on Ir(lll). By a careful study, the authors have ruled out the individual atomic 
mechanisms discussed above as well as the dislocat ion m echanism. Instead, they suggest that gliding of the cluster 
as a whole is likely to explain the observed motion [157|. 

Hamilton later studied the case of heteroepitaxial, strained islands [158]. He has shown that - due to the misfit 
between the substrate and the island structures - there can exist islands for which introducing a dislocation does not 
cost too much extra energy. These metastable mi sfit d islocations would propagate easily within the islands, leading 
to "magic" island sizes with a very high mobility [158|. 



2. 3D island diffusion mechanisms 



For the 3D clusters, the three microscopic mechanisms presented above are possible in principle. However, as 
noted above, the individual atom mechanisms lead to a diffusivity smaller than the diffusion of 5&2300 on graphite 
by several order s of magnitude. These mechanisms have also been ruled out for the diffusion of gold crystallites on 
ionic substrates |145| . Se veral te ntati ve explanations based on the gliding of the cluster as a whole over the substrate 
have been proposed [143]. Reiss |159[ showed that, for a rigid crystallite which is not in epitaxy on the substrate, the 
activation energy for rotations might be weak, simply because during a rotation, the energy needed by atoms that 
have to climb up a barrier is partially offset by the atoms going into more stable positions. Therefore, the barrier for 
island diffus ion is of the same order as that for an atom, as long as the island does not reach an epitaxial orientation. 
Kern et al. [143, 160 allowed for a partial rearrangement of the interface between the island and the substrate when 
there is a misfit. The interface would be composed of periodically disposed zones in registry with the substrate, 
surrounded with perturbed ("amorphous") zones, weakly bound to the substrate. This theory - similar in spirit to 
the dislocation theory proposed by Hamilton [ 156, 158 for 2D islands - leads to reasonable predictions [143] but is 
difficult to test quantitatively. 

To clarify the microscopi c me chanisms of 3D cluster diffusion, I now present in detail Molecular Dynamics (MD) 
studies carried out recently [ 16 1 1 . These simulations aimed at clarifying the generic aspects of the question rather than 
modeling a particular case. Both the cluster and the substrate are made up of Lennard- Jones atoms [162|, interacting 
through potentials of the form : 
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nr)=4e((f) 12 -(f) 6 ). 

Empirical potentials of this type, originally developed for the description of inert gases, are now commonly used 
to model generic properties of condensed systems. Lennard- Jones potentials include only pair atom-atom interaction 
and ensure a repulsive interaction at small atomic distances and an attractive interaction at longer distances, the 
distance scale being fixed by a and the energy scale by e. For a more detailed discussion of the different interatomic 
potentials available for MD simulations and their respective advantages and limitations, see Ref. [163]. The substrate 
is modeled by a single layer of atoms on a triangular lattice, attached to their equilibrium sites by weak harmonic 
springs that preserve surface cohesion. The Lennard- Jones parameters for cluster atoms, substrate atoms and for the 
interaction between the substrate and the cluster atoms are respectively (e cc , <J CC ), (e ss , a ss ) and (e sc , a sc ). e cc and a cc 
are used as units of energy and length. e sc , a ss and T, the temperature of the substrate, are the control parameters of 
the simulation. The last two parameters are then constructed by following the standard combination rules : e ss = a® s 
and a sc = \ (a cc + <t S s)- Finally, the unit of time is defined as r = {Mo 2 c / ecc) 1 / 2 , where M is the mass of the atoms 
which is identical for cluster and substrate atom s. The simulation uses a standard molecular dynamics technique with 
thermostatting of the surface temperature [164]. 



In these simulations, the clusters take the spherical cap shape of a soli d droplet (Fig. 29) partially wetting the 
substrate. The contact angle, which can be defined following reference [165|, is roughly independent of the cluster size 
(characterized by its number of atoms n, for 50 < n < 500. This angle can be tuned by changing the cluster-substrate 
interaction. For large enough e sc , total wetting is observed. The results presented below have been obtained at a 
reduced temperature of 0.3 for which the cluster is solid. This is clearly visible in Fig. where the upper and lower 
halves of the cluster, colored white and grey at the beginning of the run, clearly retain their identity after the cluster 
center of mass has moved over 3 lattice parameters. Hence the cluster motion appears to be controlled by collective 
motions of the cluster as a whole rather than by single atomic jumps. 

The MD simulations have confirmed that one of the most important parameters for determining the cluster diffusion 
constant is the ratio of the cluster lattice parameter to the substrate lattice parameter. The results for the diffusion 
coefficient are shown in Fig. |30|a. When the substrate and cluster are commensurate (o~ ss = a cc = 1), the cluster 
can lock into a low energy epitaxial configuration. A global translation of the cluster would imply overcoming an 
energy barrier scaling as n 2 / 3 , the contact area between the cluster and the substrate. In that case diffusion will 
be very slow, unobservable on the time scale of the MD simulations. What is interesting to note is that even small 
deviations from this commensurate case lead to a measurable diffusion on the time scale of the MD runs. This can be 
understood from the fact that the effective potential in which the center of mass moves is much weaker, as the cluster 
atoms, constrained to their lattice sites inside the rigid solid cluster, are unable to adjust to the substrate potential 
(see above, Reiss model ]159[ ). The effect is rather spectacular : a 10% change on the lattice parameter induces an 
increase of the diffusion coefficient by several orders of magnitude. 

Finally, I show in Fig. |30| b the effect of cluster size on the diffusion constant for different lattice parameter values. 
As the number n of atoms in the cluster is varied between n = 10 and n — 500, the diffusion constant decreases, 
roughly following a power law D ~ n a . This power law exponent a depends significantly on the mismatch between 
the cluster and the substrate lattice parameters. For high mismatches (a ss — 0.7,0.8), a is close to —0.66. As the 
diffusion constant is inversely proportional to the cluster-substrate friction coefficient, this result is in agreement with 
a simple "surface of contact" argument yielding D ~ n~ 2 / 3 . On the other hand, when the lattice mismatch is equal to 
0.9, one obtains a ~ —1-4, although the shape of the cluster, characterized by the contact angle, does not appreciably 
change. It is instructive to follow the trajectory followed by the cluster center of mass (Fig. |3l]). In the runs with 
a large mismatch (Fig. |3l|a), this trajectory is "brownian-like" , with no apparent influence of the substrate. This 
is consistent with the simple "surface of contact" argument. Instead, when the mismatch is small (Fig. |3l|b), the 
center of mass of the cluster follows a " hopping-like" trajectory, jumping from site to site on the honeycomb lattice 
defined by the substrate. When cr ss = \/3/2, there seems to be a transition between the two regimes around n = 200. 

It is interesting to consider the interpretation of cluster motion in terms of dislocation displacement within the 



cluster, a mech anism which has been proposed to explain rapid 2D cluster diffusion [156,158| (see the discussion in 
Section VII A 1). For this, one can "freeze" the internal degrees of freedom of the cluster deposite d on a thermalized 
substrate. The center of mass trajectory is integrated using the quater nion algorithm 161 , 164 ], Surprisingly, the 



diffusion constant follows the same power law as in the free cluster case [161]. This result proves that the diffusion 



mechanism in this case cannot be simply explained in terms of dislocation migration within the cluster as proposed 
to explain the diffusion of 2D islands in [ 156,158). As the substrate atoms are tethered to their lattice site, strong 
elastic deformations or dislocations within the lattice are also excluded. Hence, the motor for diffusion is here the 
vibrational motion of the substrate, and its efficiency appears to be comparable to that of the internal cluster modes. 

Very recently, U. Landmann performed MD simulations of diffusion of large gold clusters on HOPG substrates 
|166| . He finds high cluster mobility, in agreement with the preceding simulations. His studies show that cluster 
diffusion in this case proceeds by two different mechanisms : long (several cluster diameters) linear "flights" separated 
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by relatively slow diffusive motion as observed in the preceding simulations. Further work is needed to ascertain the 
atomic mechanisms leading to this kind of motion. 



3. Discussion 



What are the (partial) conclusions which can be drawn from these studies of cluster diffusion? I think that the 
main parameter determining the mobility of 3D islands on a substrate is the possible epitaxy of the cluster on the 
substrate. Indeed, if the island reaches an epitaxial orientation, it is likely to have a mobility limited by the individual 
atomic movements, which give a small diffusion constant (of order 10 _17 cm 2 s _1 at room temperature). Diffusivities 
of this magnitude will not affect the growth of cluster films during typical deposition times, and clusters can be 
considered immobile. The effect of these kind of diffusio n rat es can only be seen by annealing the substrates at higher 
temperatures or for long times. According to Hamilton |l56fl , dislocations could propagate even for epitaxial islands, 
but it is likely that this mechanism is more important in the case of heteroepitaxial islands which I now proceed to 
discuss. Indeed, if the island is not in epitaxy on the substrate, high mobilities can be observed because the cluster 
sees a potential profile which is not very different from that seen by a single atom. It should be noted that this non- 
epitaxy can be obtained when the two lattice parameters (of the substrate and the island) are very different, or also 
when they are compatible if there is relative misorientation. The latter has been observed for gold on ionic substrates 
[143 1 and mobility is relatively high until the crystallites reach epitaxy. The MD simulations presented above show 
that, for Lennard- Jones potentials, only homoepitaxy prevents clusters from moving rapidly on a surface. It should 
be noted that relaxation of the cluster or the substrate - which would favor a locking of the cluster in an energetically 
favorable position at the expense of some elastic energy - has not been observed in these LJ simulations, nor has 
dislocation propagation. This is probably realistic for the low interaction energies which correspond to metal clusters 
on graphite. It could also be argued that dislocation motion is more difficult in 3D clusters than in 2D islands since 
the upper part of the particle (absent in 2D islands) tends to keep a fixed structure. Another important parameter 
is the cluster-substrate interaction : one can think that a large attractive interaction (for metal on metal systems for 
example) can induce an epitaxial orientation and prevent the cluster from diffusing, even in the heteroepitaxial case. 
The differences between the diffusion of clusters grown on a substrate by atom deposition and aggregation and those 
previously formed in a beam and deposited must also be investigated. One could anticipate that islands formed by 
atom aggregation on the substrate would accommodate easily to the substrate geometry, whereas preformed clusters 
might keep their (metastable) configuration. However, it is not at all clear that island nucleation and epitaxy are 
simultaneous phenomena, for it has been observed that islands can form in a somewhat arbitrary configuration and 
subsequently orient on the substrate after diffusion and rotation (see Ref. [143]). 



B. Cluster-Cluster coalescence 

What happens now when two clusters meet? If they remain simply juxtaposed, morphologies similar to Fig. |l8| a 
are observed. In this case, the incident clusters have retained their original morphology, and the supported particles 
are identical to them, even if they are in contact with many others after cluster diffusion. It is clear, by looking for 
example to Fig. ^ that this is not always the case. In these cases, the supported islands are clearly larger than the 
incident clusters : some coalescence has taken place. How can one understand and predict the size of the supported 
particles? Which are the relevant microscopic parameters? This section tries to answer these questions, which are of 
dramatic interest for catalysis, since the activit y of the deposits crucially depends on its specific area and therefore 
on the sintering process (See for example Refs. |168| ). 

I will first briefly examine the classical theory for sphere-sphere coalescence (i.e. ignoring the effect of the substrate) 
and then review recent molecular dynamics simulations which suggest that this classical theory may not be entirely 
satisfactory for nanoparticles. 



1. Continuum theory of coalescence 



The standard analysis of kinetics of sintering is due to Mullins and Nichols [169,170]. The "motor" of the coalescence 
is the diffusion of atoms of the cluster (or island) surface from the regions of high curvature (where they have less 
neighb ors a nd therefore are less bound) towards the regions of lower curvature. The precise equation for the atom 
flux is |70j 
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•h = , „ V S A 



(7.2) 



where D s is the surface diffusion constant (supposed to be isotropic), 7 the surface energy (supposed to be isotropic 
too), Q the atomic volume, v the number of atoms per unit surface area, ks Boltzmann's constant, T the temperature 
and K the surface curvature {K=l/R\ + I/R2) where i?i and R2 are the principal radii of curva ture) . For sphere- 
sphere coalescence, an order of magnitude estimation of the shape changes induced by this flux is |170| : 



dn n d 2 K, 
-~2B^(y = S = 0) 



(7.3) 



where dn is the outward normal distance traveled by a surface element during dt, s the arc length and B 
D s^jV^v IksT (the z axis is taken as the axis of revolution). For this geometry, Eq. 7.2 becomes (Fig. |32l) : 



dl B ( I 
di~l?y~R 

where I have made an order of magnitude estimation of the second derivativ e of the curvature 
K{l))/l and similarly d 2 K/ds 2 - (1 - l/R)/l 3 (see Fig. ||). Integrating Eq. [Ml leads to 



(7.4) 



dK/ds ~ {K(R) 



I ~ (r 4 + 4Bt) 1/4 for 1<^R 



(7.5) 



Eq. [7lj gives an estimation of the coalescence kinetics for two spheres of radius r and R. 

However, despite its pla usib ility, Eq. 7J5 has to be used with care. Fir st, t he calculation leading to it from the 

assumes isotropic s urfac e tension and 
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expression of the flux (Eq. 7.2) is only approximate. More importantly, Eq 
diffusion coefficients. While this approximation may be fruitful for large particl es (i n the nm range, 171]), it is clearly 
wrong for clusters in the nanometer range. These are generally facetted [ 172| -174| as a result of anisotropic surface 
energies. This has two important consequences : first, since the particles are not spherical, the atoms do not feel a 
uniform curvature. For those located on the planar facets, the curvature is even 0, meaning that they will not tend to 
move away spontaneously. This effect should significantly reduce the atomic flux. Second, the diffusion is hampered 
by the edges between the facets [175 1 which induce a kind of " Schwoebel" effect [176|. Then, the effective mass transfer 
from one end of the cluster to the other may be significantly lower than expected from the isotropic curvatures used in 
Eq. 7.2. For these anisotropic surfaces, a more general formula which takes into account the dependence of 7 on the 
crystallographic orientation should be used (see for example Ref. However, this formula is of limited practical 

interest for two reasons. First, the precise dependence of the surface energy on the crystallographic orientation is 
difficult to obtain. Second, as a system of two touching facetted clusters does not in gene ral s how any symmetry, the 
solution to the differential equation is hard to find. One possibility currently explored [177] is to assume a simple 
analytical equation for the anisotropy of 2D islands and integrate numerically the full (anisotropic) Muffins' equations. 



2. Molecular Dynamics simulations of coalescence 



Since continuum theories face difficulties in characterizing the evolution of nan oparticle coa lescence, it might be 
useful to perform molecular dynamics (MD) studies of this problem. Several studies [172,178,179] have been performed, 
showing that two disti nct and generally subsequent pr ocesses l ead to coalescence for particles in the nanometer range 
: plastic deformation |179[ | and slow surface diffusion [ 172| , 178] . 

Zhu and Averback |17S[ have studied the first stages (up to 160 ps) of the coalescence of two single-crystal copper 
nanoparticles (diameter 4.8 nm). Fig. |3^ presents four stages of the coalescence process, demonstrating that plastic 
deformation takes place (see the arrows indicating the sliding planes) and that a relative rotation of the particles 
occurs during this plas tic d eformation (c and d). During the first 5 ps, the deformation is elastic, until the elastic 
limit (roughly 0.8 nm [179Q is reached : after this, since the shear stress (Fig. is very high, dislocations are 
formed and glide on (111) planes in the < 110 > direction, as usually seen in fee systems. Fig. [34] also shows that 
after 40ps (i.e. Fig. |33|c) the stress on the glide plane is much smaller and dislocation motion is less important : the 
two particles rotate until a low-energy grain boundary is found (Fig. 33 i). This intial stage of the coalescence, where 



the two particles reorient and find a low-energy configuration, is very rapid, but does not lead in general to thorough 
coalescence. An interesting exception might have been found by Yu and Duxbury |178| : their MD simulations 
showed that for very small clusters (typically less than 200 atoms) coalescence is abrupt provided the temperature is 
sufficiently close to the melting temperature. They argue that this is due to a (not specified) "nucleation process" : 
plastic deformation is a tempting possibility. 
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For larger clusters, the subsequent stages are much slower and imply a different mechanism : atom diffusion on 
the surf ace of the particles. The intial stages of this diffusion-mediated coalescence have been studied by Lewis et 
al. [172 1 . The point was to study if Mullins' (con tinuu m) predictions were useful in this size domain. In Lewis et al. 
simulations, the embedded-atom method (EAM) [167] was used to simulate the behavior of unsupported gold clusters 
for relatively long times (~ 10ns). Evidently, an important role of the substrate in the actual coalescence of supported 
clust ers is to ensure thermalization, which is taken care of here by coupling the system to a fictitious "thermostat" 
[ 180 1 . One therefore expects these coalescence events to be relevant to the study of supported clusters in the case 
where they are loosely bound to the substrate, e.g., gold clusters on a graphite substrate. Strong interaction of the 
clusters with the substrate may be complicated and l ead to c luster deformation even for clusters deposited at low 
energies, for example if the cluster wets the substrate [181,182]. 

Fig. [55] shows the evolution of the ratio x/R, where x is the radius of the neck between the two particles. After 
an extremely rapid approach of the two clusters due to the mechanisms studied above (plastic deformation), a slow 
relaxation to the spherical shape starts (Fig. |3^). The time scale for the slow sphericization process is difficult 
to estimate from Fig. [35] but it would appear to be of the order of a few hundred ns or more. This number is 
substantially larger th an o ne would expect on the basis of phenomenological theories of the coalescence of two soft 
spheres. Indeed, Ref. 1 1 70 1 predicts a coalescence time for two identical spheres t c = ksTR 4 /(CDgja 4 ), where D s 
is again the surface diffusion constant, a the ato mic size, 7 the surface energy, R the initial cluster radius, and C 
a numerical cons tant (C = 25 according to Ref. [ 170 1 ) ; taking D s ~ 5 10~ 10 m 2 s _1 (the average value found in the 
R = 30 A, 



172]), 
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Um 2 , and a = 3 A, this yields a coalescence time r c of the order of 40 



simulations, see 

ns. The same theories, in addition, make definite predictions on the evolution of th e sha pe of the system with time. 
In particular, in the tangent-sphere model, the evolution of the ratio x/R is found |17C] to vary as x/R ^ (t/Yc) 1 / 6 
for values of x/R smaller than the limiting value 2 1 / 3 . In Fig. [55], the prediction of this simple model (full line) is 
compared with the results of the present simulations. There is no agreement between model and simulations. The 
much longer coalescence time observed has been attributed [ 172 1 to the presence of facets on the initial clusters, which 
persist (and rearrange) during coalescence. The facets can be seen in the initial and intermediate configurations of 
the system in Fig. |36j the final configuration of Fig. 36 shows that the cluster is more spherical (at least from this 
viewpoint), and that new facets are forming. That diffusion is slow can in fact be seen from Fig. pq: even after 10 ns, 
at a temperature which is only about 200 degrees below melting for a cluster of this size, only very few atoms have 
managed to diffuse a significant distance away from the contact region. 

The precise role of the facets in the coalescence process is a subject of current interest. Experiments have shown 
that shape evolution is very slow in p resence of facets for 3D crystallites (see for example ]183| ) and recent experiments 
[184 1 and computer simulations [185] on 2D islands suggest that the presence of facets can be effective in slowing down 
the coalescence process. Clearly, more work is needed to get a quantitative understanding of nanoparticles coalescence, 
and to evaluate the usefulness of Mullins' approach, especially if one manages to include the crystalline anisotropy 
(see also Refs. [186Q. 



C. Island morphology 

Now I can turn on to the prediction of one of the essential characteristics of cluster films : the size of the supported 
particles. As I have already mentioned in the introduction, the size of the nanoparticles controls many interesting 
properties of the films. Therefore, even an approximate result may be useful, and this is what I obtain in this section. 

The experiments shown above demonstrate that the supported particles can have a variety of sizes, from that of 



the inc ident clusters (S'62300/HOPG, Section VIA) up to many times this size (for example ^4zi25o/HOPG, Section 
VI B 3). To understand how the size of the supported particles is determined, one can look at a large circular island 
to which clusters are arriving by diffusion on the substrate (Fig. [32]). There are two antagonist effects at play here. 
One is given by thermodynamics, which commands that the system should try to minimize its surface (free) energy. 
Therefore one expects the clusters touching an island to coalesce with it, leading to compact (spherical) domains. 
The other process, driving the system away from this minimization is the continuous arrival of clusters on the island 
edge. This kinetic effect tends to form ramified islands. What is the result of this competition? Since there is a 
kinctically driven ramification process, it is essent ial t o take into account the kinetics of cluster-cluster coalescence, 



as sketched in the previous section. I will use Eq. 7.5 even if it is only approximate, to derive an upper limit for the 
size of the compact domains grown by cluster deposition. It is an upper limit sin ce, a s pointed out in the previous 
section, coalescence for facetted particles could be slower than predicted by Eq. |7.5| , hence diminishing the actual 
size of these domains. 

We first need an estimate of the kinetics of the second process : the impinging of clusters on the big island. A very 



simple argument is used here (see also [187| for a similar analysis for atomic growth) : since the number of clusters 
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reaching the surface is F per unit surface per second and the total number of islands is N t per unit surface, each 
island receives in average a cluster every t r = N t /F. 

We are now in a position to quantify the degree of coalescence in a given growth experiment. Let me suppose that 
a cluster touches a large islan d at t=0. If no cluster impinges on the island before this cluster completely coalesces 
(in a time r c according to Eq. 7.5), then the islands are compact (circular). Instead, if a cluster touches the previous 



cluster before its total coalescence has taken place, it will almost freeze up the coalescence of the previous cluster. 
The reason is that now the atoms on the (formerly) outer surface of the first cluster do not feel curvature since they 
have neighbors on the second cluster. The mobile atoms are now those of the second cluster (see Fig. |3^, label a) 
and the coalescence takes a longer time to proceed (the atoms are farther form the big island). Then, if t r <C r c , 
the islands formed on t he surface are ramified. For intermediate cases, the size R c of the compact domains can be 
estimated from Eq. 7.5 as R c = x{t r ) where t r takes into account the fact that, to freeze the coalescence of a previous 
ras to touch the island at roughly the same point : t r ~ t r 2nR/r and 



cluster, one cluster 
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Eq. 7.6 describes the limiting cases {B ~ oo or B ~ 0) correctly. The problem with the interme diat e cases is to 
obtain a reliable estimate of the (average) atomic surface self diffusion. For gold, Chang and Thiel [ 136 1 give values 
which vary between 0.02 eV on compact facets and 0.8 eV on more open surfaces. One solution is to go the other way 
round and estimate D s from the experimental data and Eq. [7.6| . From Fig. ^2], estimating R c from the thickness 
of the island arms, and using the experimental values for r (0.85 nm) and the fact that since the flux is pulsed (see 
Section VI B 3), the time between two successive arrivals of clusters is approximately the time between two pulses (0.1 
s), and not N t /F, one obtains D s ~ 3 10~ 3 cm 2 s~ x exp(— 0.69eV/ (ksT)), which seems a sensible value. 

Despite the difficulty of defining average diffusion coefficients, one can use Eq. |7.6j to obtain a reasonable guess 
for the size of the compact domains by assuming that D s is thermally activated : D S (T) = Doexp(— E a / (ksT)) with 
a prefactor Do = 10 _3 cm 2 s _1 and an activation energy E a taken as a fraction of the bonding energy between atoms 
(proportional to &bT/). One obtains [18S] 



B = 10 11 exp(-4.6T / /r)nm 4 /s 



(7.7) 



Inserting this value in Eq. 7.£ leads to Fig. [57] where the size of the compact domains is plotted as a function of 
T/Tf. The important feature is that as long as T/Tf < 1/4, the incident particles do not merge. Note that this 1/4 
is sensitive to the assumed value of D* , but only via its loga rithm. Again, this estimation of R c /r is an upper limit 
since coalescence could be slower than predicted by Eq. 7.5. 



What would happen now if the incident clusters were liquid! An experimental example of this liquid coalescence is 
given by t he de position of Tnioo on a-C (see above). A rough guess of the coalescence time is given by a hydrodynamics 
argument [189|: the driving force of the deformation is the surface curvature j/R 2 where 7 is the liquid surface tension 
and R the cluster radius. This creates a velocity field which one can estimate using the Navier-Stokes equation : 
r]Av — 7/i? 2 where 77 is the viscosity and v is the velocity of the fluid. This leads to T c (liquid) ~ R/v ~ r]R/j. Inserting 
reasonable values for both rj (.01 Pa s) and 7 (1 J mT 2 ) leads to r c (liquid) ~ O.Oli? which gives T c (liquid) = lOps for 
R ~ Inm. This is t he g ood order of magnitude of the coalescence times found in simulations of liquid gold clusters 
(T c (liquid) ~ 80ps |l72[ |). Now, since T c (liquid) -C t r (t r ~ 0.1s, see above) cluster-cluster coalescence is almost 
instantaneous, which would lead to R c ~ 00. In fact, R c is limited in this case by static coalescence between the big 
islands formed during the growth. The reason is that the big islands may be solid or pinne d by defects leading to a 



slow coalescence. The analysis is similar here to what has been done for atomic deposition [19C] 



D. Thick films 



The preceding section has studied the first stages of the growth, the submonolayer regime, which interests researchers 
trying to build nanostructures on the surface. I attempt here a shorter study of the growth of thick films, which 
are known to be very different from the bulk material in some cases [^,^,^3|. The main reason for this is their 
nanostructuration, as a random stacking of nanometer size crystallites. Therefore, it is interesting to understand 
how the size of these crystallites is determined and how stable the nanostructured film is. One can anticipate that 
the physical mechanism for cluster size evolution is, as in the submonolayer case, sintering by atomic diffusion. For 
thick films however, surface diffusion can only be effective before a given cluster has been "buried" by the subsequent 
deposited clusters. Thus, most of the size evolution takes place during growth, for after the physical routes to 
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co alesc ence (bulk or grain boundary diffusion) are expected to be much slower. Studies of compacted nanopowders 
|^,192 1 have shown that nanoparticles are very stable against grain growth. Siegel explains this phenomenon 



in the following way. The two factors affecting the chemical potential of the atoms, and potentially leading to 
structure evolutions are local differences in cluster size or in curvature. However, for the relatively uniform grain size 
distributions and flat grain boundaries observed for cluster assembled materials ||, these two factors are not active, 
and there is nothing telling locally to the atoms in which direction to migrate to reduce the global energy. Therefore, 
the whole structure is likely to be in a deep local (metastable) minimum in energy, as observe d in closed-cell foams. 



The stability of such structures has been confirmed by several computer simulations [193,194] which have indicated 



a possibl e m echanism of grain growth at very high (T/Tf ~ 0.8) temperatures : grain boundary amorphization or 



melting [ 193 1 



What determines the size of the supported particles during the growth? For thick films, a reasonable assumption 
is that a cluster impinging on a surface already covered by layer of clusters does not diffuse, because it forms strong 
bonds with the layer of the deposited clusters. This hypothesis has been checked for the growth of 562300 on graphite 



[ 1 1 5 1 . There are two main differences with submonolayer growth : first, an impinging cluster has more than one 
neighbor and the sphere-sphere kinetics is not very realistic; second the relevant time for ramification I have used in 
the preceding section is no more useful here since clusters do not move. As a first approximation, to obtain an upper 
limit in the size of the domains, we can use the same coalescence kinetics and take a different "ramification" time : 
the average time for the arrival of a cl uste r touching another is roughly tf ~ \/(Fd 2 ) where d — 2r is the diameter of 
the cluster. If the same formula (Eq. ff^ ) is used, one finds 

*c 4 = r 4 + ^ (7.8) 
The results obtained using the same approximation as in the preceding section for B (Eq. |7.7|) are shown in Fig. 



Experimentally, there are observations for deposition of Ni and Co clusters [ 191 [. The size of the crystallites is 
comparable to the size of the incident (free ) cl usters. This is compatible with Eq. 7.S since the Tf of these elements 
is very high (~ 1800i<Q. Therefore, Eq. [T| predicts that films grown at T = 300K, (T/T f ~ 0.17) should keep 



a nanostructuration with R c ~ r, as is observed experimentally [191]. I stress again that a structure obtained with 
cluster deposition with this characteristic size is not likely to recr ystallize in the bulk phase (thereby loosing its 
nanophase properties) unless brought to temperatures close to Tf p|, 193| ,194]. 



VIII. CONCLUSIONS, PERSPECTIVES 

What are the principal ideas presented in this paper? 

First, useful models to analyze the first stages of thin film growth by cluster deposition have been presented in 



detail (Section |nj). These models are useful at a fundamental level, and I have shown in Section VI how many 
experimental results concerning submonolayer growth can be interpreted by combining these few simple processes 
(deposition, diffusion, evaporation . . . ). Specifically, by comparing the experimental evolution of the island density as 
a function of the number of deposited particles to the predictions of computer simulations, one can obtain quantitative 
information about the relevant elementary processes. 

Second, the quantitative information on diffusion has shown that large clusters can move rapidly on the surface, 
with diffusion constants compa rabl e to the atomic ones. A first attempt to understand this high diffusivity at the 



atomic level is given in Section VII : the conclusion is that rapid cluster diffusion might be quite common, provided 



the cluster and the substrate do not find an epitaxial arrangement. Concerning cluster-cluster coalescence, it has been 



suggested that this process can be much slower than predicted by the usual sintering theories |170|, probably because 
of the cluster facets. 

Third, despite all the approximations involved in its derivation, Fig. ^ gives an important information on the 
morphology of the film : an upper limit for the ratio of the size of the compact domains over the size of the incident 
clusters. This helps understanding why cluster deposition leads to nanostructured films provided the deposition 
temperature is low compared to the fusion temperature of the material deposited (T s < Tf /4) . Clearly, further 
experimental and theoretical work is needed in order to confirm (or invalidate) Fig. B7l 

Ii is clear that we still need to understand many aspects of the physics of cluster deposition. Possible investigation 
directions include the following, given in an arbitrary order. First, the coalescence of nanoparticles has yet to be 
understood and quantified. This is a basic question for both submonolayer and thick materials. Second, one has 
to characterize better the interactions between clusters and the substrate, and especially its influence on the cluster 
diffusion. It is also important to investigate the possible interactions between the clusters, which could dramatically 
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affect the growth. Obtaining ordered arrays of nanoparticles is a hot topic at this moment. A possibility is the 
pinning of the clusters on surface " defects" which demands a better understanding of cluster interaction with them. 
Another idea is to use t he s elf-organization of some bacteria to produce an ordered array on which one could arrange 



the clusters (See Ref. [ 195 1 , especially Chap. 5). Clearly investigating the interaction of clusters with biological 



substrates is not an easy task, but it is known that practical results are not always linked to a clear understanding of 
the underlying mechanisms . . . 
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FIG. 1. Ag nanoislands grown on two monolayers of Ag deposited on Pt(lll) and annealed at 800K. The inset shows a fast 
fourier transform of the spatial distribution. From Ref. /mj/ ) 

FIG. 2. By changing the mean size of the incident antimony clusters, one can dramatically change the morphology of the 
submonolayer film. The four micrographs have been obtained for the same thickness (1 nm) and deposition rate (5x 10~ 3 nms~ 1 ). 
The mean sizes are : (a) Sb^, (b) Sbi6, (a) Sb^G, (a) Si^O. The changes in morphology are interpreted by the different 
mobilities of the clusters as a function of their size, as well as their different coalescence dynamics and sensitivity to surface 
defects. From Ref. JIj[ /. 

FIG. 3. Molecular Dynamics simulations of the morphology of films obtained by M01043 cluster dep osition for increasing 
incident kinetic energies per atom (as indicated in the figures) onto a Mo(001) substrate. From Ref. pqj 
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FIG. 4. Principle of a Kinetic Monte Carlo simulation (see text). 



FIG. 5. Mam elementary processes considered in this paper for the growth of films by cluster deposition, (a) adsorption of a 
cluster by deposition; (b) and (d) diffusion of the isolated clusters on the substrate; (c) formation of an island of two monomers 
by juxtaposition of two monomers (nucleation) (d) growth of a supported island by incorporation of a diffusing cluster (e) 
evaporation of an adsorbed cluster. I also briefly consider the influence of island diffusion (f). 



FIG. 6. Possible interaction of two clusters touching on the surface : (a) pure juxtaposition (b) total coalescence. Interme- 
diate cases (partial coalescence) are possible and will be described later. 



FIG. 7. Time scales of some elementary processes considered in this paper for the growth of films by cluster deposition. The 
relevant processes are those whose time scale are smaller than the deposition time scale shown by the arrow in the left. In this 
case, models including only cluster diffusion on the substrate and cluster cluster coalescence are appropriate. "Island diffusion" 
refers to the motion of islands of clusters as a whole, "cluster dissociation" to the evaporation of atoms from the cluster and 
"interdiffusion" to the exchange of atoms in the cluster with substrate atoms. 



FIG. 8. Evolution of the monomer and island densities as a function of the thickness (in monolayers), for islands formed 
by pure juxtaposition : (a) complete condensation, F = 1CP 8 , T e = 10 10 (r = 1). These values mean : Xs = 10 5 and ice = 22 
(b) important evaporation, F = 10 -8 , r e = 600 (r = 1) (X s = 2 5 and ice = 11). £cc represents the mean island separation at 



saturation for the given fluxes when there is no evaporation (117]. The length units correspond to the incident cluster (monomer) 



diameter. In (b) the "condensation" curve represents the total number of particles actually present on the surface divided by the 
total number of particles sent on the surface (F t). It would be 1 for the complete condensation case, neglecting the monomers 
that are deposited on top of the islands. The solid line represents the constant value expected for the monomer concentration, 
while the dashed line corresponds to the linear increase of the island density (see text). 



FIG. 9. Morphology of a submonolayer deposit in the case of growth with complete condensation and pure juxtaposition : 
(a) 9 = 0.1% ; (b) 9 = 15%. The values of the parameters are : F = 10~ 8 ML/s, r = 1 and L=300. 



FIG. 10. Saturation island density as a function of the normalized flux (r = 1) for different growth hypothesis indicated 
on the figure, always in the case of island growth by pure juxtaposition, "no evap" (circles) means complete condensation. 
Triangles show the densities obtained if there is evaporation, for r e = 100. In the preceding cases, islands are supposed to be 
immobile. This hypothesis is relax ed f or the last set of data, mobile islands (squares) , where island mobility is supposed to 
decrease as the inverse island size [102 (there is no evaporation). The dashed line is an extrapolation of the data for the low 
normalized fluxes. Fits of the different curves in the low-flux region give : "no evap" (solid line): N sat = 0.53(_Ft) 0,36 ; "evap" 
(dotted line) : N 3a t = 0.26-F°' 67 r _1 ^ 3 r e (for the r and r e exponents, see 117| and Appendix A) and "mobile islands" (dashed 
line) : N sat = 0.33(Fr) 042 



FIG. 11. Evolution of the island density as a function of the thickness for different growth hypothesis. This figure shows that 
the same saturation density can be obtained for films grown in very different conditions. The different sets of data represent 
: triangles : growth with coalescence and evaporation, r e = lOOr and Ft = 1.2 10 -8 , circles : growth with coalescence but 
without evaporation (Ft = 3 10 -10 ), solid line : growth with pure juxtaposition without evaporation (Ft = 2.5 10~ 9 ), squares 
: growth with coalescence on defects (defect concentration : 5 10 -4 per site) and Ft = 10~ 14 (no evaporation), dashed line : 
growth with pure juxtaposition without evaporation but with mobile islands, Ft = 10~ 8 . 



FIG. 12. Normalized island size distributions for F = 10 , r = 1 and different values of r e for islands formed by pure 
juxtaposition (no coalescence). The size distributions were obtained for different coverages 9 between .05 and 0.2. The solid 
line shows the size distribution obtained without evaporation, the dashed line that obtained with mobile dimers and the 
numbers show the different values of T e . The size distributions shown here have been obtained with Ft = 10~ 8 , but the same 
distr ibutions are obtained if both F, r and r e are changed but the same parameter e = (1 + Xs)X s (Ft) is obtained (see Ref. 



30 



FIG. 13. Saturation island density as a function of the normalized flux (r = 1) for different growth hypothesis in the case of 
growth by total coalescence (3d islands). I show the densities obtained for the complete condensation case (filled circles) and 
for two different evaporation times : r e = 100 (triangles) and r e = 20 (squares). The label defects means growth in presence of 
defects which act as nucleation centers. Their concentration is 10 -3 per site. The dashed line is an extrapolation of the defect 
data for the low normalized fluxes. Fitting the simulation data leads to the following numerical relations : N sat = 0.27(Ft)°' 2S6 
when there is no evaporation (solid line) ; N sa t ~ 0.039i ?0 ' 55 T~ 2//3 Te when evaporation is significant (from an approximation 
for the two dotted curves) : the exponents for r and r e have been derived from a rate-equations treatment (Appendix B) 

FIG. 14. Normalized island size distributions obtained for F — 10 -8 , r = 1 and different values of the evaporation time r e 
for islands formed by total coalescence. The size distributions were obtained for different coverages 9 between .05 and 0.2. The 
solid line shows the size distribution obtained without evaporation. The number next to each symbol corresponds to r e . 

FIG. 15. Effect of the presence of defects on the island size distribution. The rescaled island size distributions are obtained 
for F — 10 -8 and different values of the evaporation time r e (r = 1) for islands formed by total coalescence by nucleation on 
defects. The size distributions were obtained for different coverages 9 between .05 and 0.15. Contrary to what is observed for 
homogeneous nucleation, (i.e. without defects) the histograms do depend on the coverage for nucleation on defects. The solid 
line shows the size distribution obtained without evaporation. 

FIG. 16. Values of (a) the thickness e aa t and (b) island density N aa t at the saturation of island density as a function of 
the e vaporation parameter n = FtX% for growth with pure juxtaposition [1111. The solid lines represent theoretical predictions 



till]. 



FIG. 17. Values of the thickness e sat (a) and the condensation coefficient C sa t (b) at the saturation of island density in the 
total coalescence limi t (3d growth for atomic deposition). In the limit of low island densities (i.e. high evaporation rates), C aa t is 



a constant (see Ref. [1341, this regime is indicated by the solid line). However, there are crossover regimes which depend on the 
precise r e and which are shown here. Then, from a measure of C sa t and N sa t one can get an estimate for r e for the not too low 
island densities which correspond to many experimental cases. In the same spirit, (a) shows the evolution of e sa t as a function 
of N a at in the crossover regime. The numbers correspond to the different t 3 /t used for the simulations and CC refers to the 
case of complete condensation (no evaporation). The dotted line in the higher left shows the limiting regime e sa t ~ Nsat -1 ■ 

FIG. 18. Typical island morphologies obtained experimentally by TEM (a) and from the computer simulations (b) at the same 
coverage, (a) S&2300 deposition on graphite HOPG at T B = 353K and f — 610~ 3 nms~ 1 , corresponding to F — 1.710 _3 ML/s. 
The deposited thickness is 0.5 nm or e — 0.14ML (b) model including only deposition, diffusion and pure juxtaposition of the 
incident clusters, Ft — 9 10 -11 

FIG. 19. (a) Evolution of the island density as a function of the deposited thickness. The solid line is a fit to the experimental 
data with Ft — 1.75 10~ 8 . (b) Evolution of the maximum island density (N sat ) as a function of the incident flux F at room 
temperature. The solid line is a fit to the experimental data : we find N sa t = a f °> 37±0 - 03 ( c ) Dependence of the diffusion 
coefficient on the temperature. From a fit on the experimental data (solid line), one finds D = Do exp (-E a /kT), with E a — 
0.7±0.1 eV and Do = 10 4 cm 2 s _1 . The island densities are expressed per site, a site occupying the projected surface of a 
cluster, equivalent to 2.08 10 _13 cm 2 

FIG. 20. Scanning electron microscopy of a submonolayer deposit of Agieo slightly accelerated (50eV) clusters on HOPG. 
From Ref. Mi. 



FIG. 21. (a) Morphology of a Sb^ film at e = 1.8ML. (b) Evolution of the island density (per site) as a function of 
thickness (ML). The dashed line represents a fit of the data with Ft — 10 -5 assuming a pyramidal (half-sphere) shape for the 
supported islands, while the solid line assumes that islands are spherical and Ft = 3 10~ 6 

FIG. 22. Morphologies of a AU250 films at e = 0.12 ML and increasing temperatures as indicated in the micrographs. There 
are less and less islands as the substrate temperature is raised and the islands become more and more compact. 
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FIG. 23. (a) Saturation island density as a function of the diffusion time t (Fi=6 ML/s) for two hypothesis : only the 
monomers move (solid line) or islands up to the pentamer move too (dashed line). The lowest island densities have been 
extrapolated, (b) Temperature dependence of diffusion coefficient as derived from (a) and Fig. ^^in the two hypothesis : only 
the monomers move (solid line) or islands up to the pentamer move too (dashed line). 

FIG. 24. Morphology (a) and diffraction pattern (b) of a AU250 submonolayer deposit on NaCl at e — 0.12 ML. The 
supported islands are small (mean diameter ~ 5nm) and in epitaxy with the substrate as shown by the diffraction pattern. 

FIG. 25. Radial distribution functions for gold clusters grown at 293K by atomic deposition before and after annealing at 
390 K for several minutes. The right inset shows the size distribution of the clust ers, which does not change, demonstrating 



that no particle-particle coalescence via atomic evaporation takes place. From Ref. t 13F. 



FIG. 26. Individual atomic mechanisms leading to island diffusion. PD refers to diffusion of atoms on the periphery of the 
island, while the exchange of atoms between the island and the atomic 2D gas is shown by the condensation and evaporation 
labels. The dashed circles represent the old positions of the atoms, while the continuous circles represent the new positions, after 
the elementary process. 

FIG. 27. Principle of island motion by dislocation propaga tion. The atomic column in the middle of the island jumps from 



fee to hep sites, moving the island center of mass. From Ref. [156] 



FIG. 28. Successive positions of a Irig 2D cluster on a Ir(lll) surface observed by field ion microscopy at low temperature. 
The motion takes place at T ~ 690if and the figures correspond to 6,10 and 14 heating intervals of 10 seconds each. From Ref. 



FIG. 29. Configuration of the Lennard- Jones cluster on the crystalline surface. (A) Side view : The cluster is partially 
wetting the surface. (B) Top view : The two halves of the cluster have been colored at the beginning of the run. After the cluster- 
center of mass has moved by roughly three substrate lattice constants from its original position, the two parts of the cluster 
are still well distinct. Then, the cluster diffusion cannot been explained in terms of single atom mechanisms (n=100, a sa =0.7, 
e ac =0.4, T=0.3) 

FIG. 30. (a) Dependence of the diffusion coefficient on the mismatch between the lattice parameter of the substrate and the 
cluster. A small change in the lattice parameter of the cluster leads to a huge change in the diffusivity (n=100, e sc =0.4, T=0.3, 
Run Length = 12500 t) (b) Dependence of the diffusion coefficient of a cluster as a function of its number of atoms. Data 
correspond to different mismatches between the cluster and the substrate lattice parameters. The diffusion coefficient decreases 
as a power law with exponent a. The two different slopes correspond to different diffusion regimes : the weaker dependence 
corresponds to a brownian trajectory; the stronger correspond to a "hopping-like" diffusion. For comparison, the arrow indicates 
the diffusion coefficient of a single adatom with a sa = 0.9. 

FIG. 31. Trajectory of a cluster center of mass diffusing on a substrate. The solid line represents the trajectory and the 
circles the equilibrium position of the surface atoms, (a) large mismatch : the motion is "brownian-like" , i.e. the cluster does 
not "see" the structure of the surface. The values of the parameters are : e sc — 0.4, a ss = 0.7, T—0.3, n=100; (b) small 
mismatch : the cluster center of mass jumps from one hexagon center to a nearest neighbor one. The values of the parameters 
are the same as for (a) except for a ss = 0.9 

FIG. 32. Schematic illustration of the competition between coalescence and kinetic ramification. R is the radius of the largest 
island, r that of the incident clusters and I stands for the typical length of a coalescing cluster. The label "a" refers to the 
ramification process when a cluster is touched by another one before coalescence can take place. 
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FIG. 33. Atomic positions at four different times during partial coalescence of Cu clusters (4-8 nm diameter each). The 
atomic positions are projected onto the (121) plane of the bottom sphere for (a) to (c) and onto the (101) plane for (d). The 
arrows indicate the sliding plane ((a) to (c)) and the grain-boundary dislocation (d). From Ref. l 17i]. 



FIG. 34. Distribution of shear stress in the sliding plane at different times during the coalescence shown in Fig. Filled 
circles : 5ps, open squares : lOps, open triangles : 20ps, crosses : SOps and plus : 40ps. From Ref. [17&I- 



FIG. 35. Evolution in ti me of the ratio of the neck radius, x, to the cluster radius, R. The full line repres ents the numerical 
solution obtained by Nichols t 17Cj with an arbitrary time scale, while the crosses are the results of Lewis et al. [17ij simulations. 



FIG. 36. Successive cluster morphologies during the coalescence of a gold 767-atom liquid cluster with another gold 1505-atom 
solid cluster. The figures represent three stages of the coalescence process after 0, 1 and 10 ns, i.e. times much longer than 
those studied in Ref. (17iJ. 



FIG. 37. Approximate dependence of the radius of the supporte d pa rticles R c as a function of the substrate temperature for 
submonolayer and thick films. Lines refer to predictions from Eq. 7. 6 with different incident cluster radius while symbols rep- 
resent experime nta l results shown. The theoretical predictions for the submonolayer regime were obt aine d by taking Nt/F — 0.1 
and using Eq. \7. \ For the thick film limit, I have taken F — 10~' A ML/s and r — 2.5nm (Eq. 7.8). One should consider 
these theoretical R c values as an upper limit since coalescence may be much slower at these (nano)scales (see the text). As a 
consequence, it is no surprise that the predicted values are clearly larger than the experimental ones. Concerning Sb, the huge 
difference can come from a partial oxidation of the clusters on the substrate because of the relatively bad vacuum conditions 
(pressure ~ 10~ 7 Torr). Even a thin oxide layer can decrease significantly atomic surface diffusion and transport, thus slowing 
the coalescence process. 
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Appendix A 

Regimes predicted by rate-equations calculations for the growth of 2D islands with evaporation. These predictions 
agree with the computer simulations presented in t his p aper and are relevant for both cluster and atomic deposition 

(see Ref. [117| for more details). 



I will here rapidly recall how the rate-equations can be written [117|, and then turn on to the different regimes 
which can be derived from them. 

The rate equation describing the time evolution of the density p of monomers on the surface is, to lowest relevant 
orders in F: 

^.=F(l-9)--^-Fp-2a oP -a i N t (8.1) 
at T e 

The first term on the right hand size denotes the flux of monomers onto the island free surface, (9 is the island coverage 
discussed below). The second term represents the effect of evaporation, i.e. monomers evaporate after an average 
time r e . The third term is due to the possibility of losing monomers by effect of direct impingement of a deposited 
monomers right beside a monomer still on the surface to form an island. This "direct impingement" term is usually 
negligible, and indeed will turn out to be very small in this particular equation, but the effect of direct impingement 
plays a crucial role in the kinetics of the system in the high evaporation regimes. The last two terms represent the 
loss of monomers by aggregation with other monomers and with isla nds respectively. The factors 0~ o and Oi are the 



"cross sections" for encounters and are calculated in Refs. [90,36 117 
The number N t of islands will be given by: 

-^±=Fp + a oP (8.2) 

where the first term represents the formation of islands due to direct impingement of deposited monomers next to 
monomers already on the surface, and the second term accounts for the formation of islands by the encounter of 
monomers diffusing on the surface. 

For the island coverage 9 i.e. the area covered by all the islands per unit area, one has: 

^- = 2 [Fp + <7 p] + <J>N t + JN t (8.3) 

The term in brackets represents the increase of coverage due to formation of islands of size 2 (i.e. formed by two 
monomers) either by direct impingement or by monomer-monomer aggregation. The next term gives the increase of 
coverage due to the growth of the islands as a result of monomers aggregating onto them by diffusion, and the last 
term represents the growth of the islands due to direct impingement of deposited monomers onto their boundary, 
or direc tly o n the island. This last term depends on Xg, the desorption length of monomers diffusing on top of the 



islands [117|. In all the simulations presented here (Section III), I have taken Xg = 0. The total surface coverage is 
given by 9 + p ~ 9 except at very short times. 

The cross sections can be evaluated in the quasistatic approximation, which consists in assuming that R does not 



vary in time and that the system is at a steady state. One finds |9(],|9(], 1 1 7 ] 



dP\ _ n / R \ K X {R/X S ) 

^)^ R = 2nDp \x- s )K^RlXs 



= 2 * RD -r = ^Dp — * ' ° (8.4) 



The cross section for monomer- monomer encounters 0~ o is obtained from the same formula substituting R by the 
monomer radius, and D by 2D as corresponds to r elati ve diffusion. 



After some additional approximations, one finds [ 117 three principal regimes which are spanned as the evaporation 
time r e decreases. They have been called : complete condensation regime where evaporation is not important, 
diffusion regime where islands grow mainly by diffusive capture of monomers and finally direct impingement regime 
where evaporation is so important that islands can grow only by capturing monomers directly from the vapor. Within 
each of these regimes, there are several subregimes characterized by the value of X s . I use lec = (^ T ) _1 ^ 6 ! the 
island-island distance at saturation when there is no evaporation and R sa t as the maximum island radius, reached at 
the onset of coalescence. 
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complete condensation X$ S> Ice 



N sat ~ J F 1 /3 T i/3 for any x * 



(8.5) 



diffusive growth 1 <C Xs <C Ice 

{ {FXlT e f/\X s +X* s )-^ if X*«i? sat (a) 

if X*»i? sat (b) 



Fr e X 2 s 



(8.6) 



withi? sat ~ (Xs+X^^iFXjre)- 1 / 3 , which gives for the crossover between regimes (a) and (b) : Xg (crossover) 
(FXlr e )-^. 



direct impingement growth Xs <C 1 



N sat ~ { 



if (a) 
(FT e ) 2 / 3 X*- 2 / 3 if l«X*«i? sat (b) 
Fr e if X* s ^R sat (c) 



5.7) 



with R sat ~ (Fre)- 1 / 3 ^ 173 , which gives for the crossover between regimes (a) and (b) : X s (crossover) 
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Appendix B 



I present here the summary of the different limits of growth of 3d islands in presence of evaporation and/or defects. 



These results are derived in detail in Ref. [ 134 1 from the resolution of rate-equations similar to those presented in 
Appendix A. For each regime, I give in the order the saturation island density A sat , the thickness at saturation e sat 
(i.e. the thickness when the island density first reaches its saturation value), the thickness at coalescence e c (i.e. the 
thickness when the island density starts to decrease due to island-island coalescence), and the scaling kinetics of the 
mean island radius as a function of time before the saturation island density is reached. I use lec = {Ft) 1 /" 7 for 3d 



islands [ 134 1 and X s = y/r e /r. 



Clean substrate (no defects) 

high evaporation : Xs <C lec "C Idef 
N sat ~ [Fr e (l + A 2 )] 2 / 3 
e sa t~e c ~ [Fr e (l + X s 2 )]-V3 

i?~ Ft 

low evaporation : l C c < As <C Idef or l C c < Idef < A s 

Nsat ~ (5) 2/? 

/mi/7 

&sat ^ € c ^ I p J 

R ~ (FDt 2 ) 1 / 9 ~ i 2 / 9 

Dirty substrate (many defects) 

high evaporation : Xs <C Idef "C lec 

N S at ~ C 

p 1 

«C c l/2 

low evaporation : l de f <C A s <C Zee or Z de / <C Zee -C X s 

N S at ~ C 

e sa t ~ c 
e c c i/ 2 

i? ~ f — ) for i < c/F, i.e. before saturation 

i? ~ (Ft/c) 1 / 3 between saturation and coalescence (c/F <t < l/Fc 1 ' 2 ). 
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Table I 

Principal symbols and terms used in this paper. The natural length unit in the model corresponds to the mean 

diameter of an incident cluster. 



Sihti hniQ nun "tcvnm e 
O y 1 I lUUlb LtllLl Ltl lllb 


J Tqn ii~ q f? ptyi nrlcc 
U IllLb, lltlllLlI tvb 


i q In n n 
lb till til 


Qftyi) //■>"/■') ; Tr\Y' r m£>ri ni~i ~tri£> a irin f£> r\ii n nnr&nnTi /in nr /'/i/e'r^ne 
bLI LLtLLil t JLfl 1 1 ItLL UIL LI It bill J lltt U if LlLftJI tUULlU 1 I UJ tlllb Lt 1 b 


71 


niiTnhpT nf ntnTYiQ nf fnp pin^fpT 

1 LLtl 1 tUC- 1 L/J UiLL/llbJ L/J 01 tC Ufc Ltd OC 1 


u. 


r> 1 i i o~l c if rlin inn cr* i n n -m rl rl _ it / ■ inn & nrc r\ ri cn&nn ri a r\nr> ~lric do inn onn ~h 

ClLlbLtl (1 Oil lie, Lt 1 ill lllll, Ll — LlQll ' LulltIC LlQ atptliab OIL Lilt, tlClllCllL 


site 


rt <r*£>n nr^r* 1 ) i otiq orl t rt r* 1 1 i ch onn nm ~lJio at T'frt oo a ~¥ o — tt rl I /I 
ill til Utt LljJ ltd U y LL tlLLbLtl Ull Lilt bill J Lltt blLt — /ILL 1 ^± 


Ml, 


tn n T) n 1 ni i r> t* • ~trir> /rm/innf nr ty~i nil or moorlorl t n r'niior i iti i inrini 111 Trip 
1 1 lUllULLiytl . Lilt LLIIIU LLIIL UJ lllllLLtl 1 ttttltLL LU tU lit I LlllljUl llliy Lilt 




bLLObLldlt WlLll Ollt LUyti UJ ClUbLtl b (J ClUbLtl [Jtl blltj 


r 


1 flip Illy Illy J I LlJu tJL,pi tbbttl III lllu 1 lULLiytl b [Ul tlLLbLtl b ptl b ILt l ptl b ttUILlh 


T 


ULJjLlblUll Lllflt . Illtilll LLIIlt lltttltU, Jul LL tlLLbLtl LU llltllvt (1 JLLIlip 




np~ti n(->(->n ftnn af(->c (in cc/Tinne 1 
UtLWttl 1 LUJU btitb 1 III ottU 1 ILlb 1 


T e 


pi^nnnrn finn firn p * rn r> n <n i~i rn p> hp inrp n td piti nrn pr pi}n r nr\ r rn~l p c irn m ~tnp qoi rin pp 
t UtipU 1 LILlUll ill I It . IlltUlt LIU It Lit J Ul t Ll IllUllUllltl t ULipUl LILto J 1 (Jill Lilt b LLI JLltt 


A S 


TfitoM diffusiun length on the substrate befuve desuvptiun : X$ ~ yj Dt c 


(p 


N uTTfiolized flux (cf) — Ft ) expressed in clusters per site 


n 
u 


7~)<7 f-p/i i on /-i i*i r-r\o fT! /~"i enr cvniT'o c o orl tin fnm c ( D — Ofho / J t ) 
ULJjLlblUll tUtJJltltllL tXpitbbtLl ill till b ( U — blLt/Jflj 


e 


mean thickness uf the film, e = Ft where t is the depusitiun time 


6 


coverage; fraction of the substrate covered by the clusters 
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the substrate over the thickness) at saturation 
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the island-island distance at saturation when there is no evaporation 
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